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We have studied the structural and electronic phase diagrams of CeFeAsOi-iF^r and 
SmFeAsOi-ajFj; by a detailed analysis of muon spin relaxation experiments, synchrotron X-ray 
diffraction, Mossbauer spectroscopy, electrical resistivity, specific heat, and magnetic susceptibility 
measurements. In these systems structural and magnetic degrees of freedom are strongly coupled 
and magnetism competes with superconductivity for the ground state. The bulk long range mag- 
netic phase transition is always preceded by a tetragonal to orthorhombic structural phase transition 
that are both suppressed on doping. However, in the Sm system we find small volumes that show 
long range magnetic order already «10 K above the structural transition temperature. In both 
systems the structural transition is absent in superconducting materials. We find a mixed phase of 
superconductivity and short range magnetic order for one composition in each system. In view of 
available literature we conclude that it remains unclear whether this is electronic phase separation on 
a nanometer length scale or microscopic coexistence. Our muon spin relaxation experiments reveal 
short range magnetic order in a broad temperature range for all magnetic materials. For the Ce sys- 
tem, short range order develops into long range order close to a percolation of magnetically ordered 
clusters. Whereas for the Sm system, the magnetically ordered clusters already show coherent muon 
spin precession, an indication for long range order. In both systems also the structural transition is 
broad, i.e., it is preceded by an increase of the Bragg peak width over a broad temperature range. 
We quantitatively analyze the broadening of the magnetic transition. In the superconducting state, 
the temperature dependence of the magnetic penetration depth is, in all studied cases, compatible 
with a single nodeless s-wave gap. We compare the magnetic and structural transition temperatures 
determined from direct {fiSH and synchrotron X-ray diffraction) and indirect (electrical resistivity) 
measurements. 

PACS numbers: 74.70.Xa, 61.05.C-, 75.30.Fv, 76.75.+i, 76.80.+y 



INTRODUCTION 

Two intensely studied phenomena in correlated elec- 
tron systems arc superconductivity (SC) and quantum 
criticality (QC). The latter arises from the suppression 
of a continuous phase transition to zero temperature at a 
quantum critical point (QCP). In its vicinity critical fluc- 
tuations play an important role for the electronic prop- 
erties of solids. Superconductivity is often found in the 
vicinity of a magnetic quantum critical point, i.e., near 
the suppression of a magnetic phase, e.g., by doping or 
pressure. It is viridely believed, that magnetic fluctua- 
tions near a magnetic QCP may play a crucial role in the 
pairing mechanism of Cooper pairs in correlated electron 
systems. Both in heavy fermion superconductors and the 



copper oxide high- Tc superconductors this interplay is 
studied. CePd2Si2[l] is an example for a heavy-Fermion 
system that is antiferromagnetic at ambient pressure, and 
becomes superconducting under pressure near the QCP 
of the antiferromagnetic order. In the cuprate supercon- 
ductor La2-2;Srj;Cu04_y quantum criticality arises from 
the suppression of the antiferromagnetic order by hole 
doping via substitution of La by Sr (for a review see 
Ref. 2 and references therein). The emergence of low 
energy spin fluctuations are believed to be crucial not 
only for the pairing mechanism of superconductivity but 
also for the electronic properties over a wide range of 
its phase diagram[3]. This general behavior under hole 
doping is found in other cuprate materials as well. 

After the discovery of the ferropnietide superconduc- 
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tor LaFeAsOi-ajFj; a plethora of different phase dia- 
grams were investigated. Among the first were the F 
doped group of rare earth i?=La, Ce, and Sm contain- 
ing i?FeAs0i_2:Fa; [4-6]. The undoped systems are semi 
metals and show a tetragonal to orthorhombic structural 
phase transition on lowering the temperature followed by 
a transition to a spin density wave phase (SDW). Upon 
doping both phases are suppressed and superconductiv- 
ity emerges with a maximal Tc~56 K. 

The phase transition from magnetic order to super- 
conductivity as a function of doping for iZ=La is abrupt, 
with only a small decrease of the magnetic ordering tem- 
perature before superconductivity emerges and magnetic 
order is fully destroyed [6]. The available experimental re- 
ports for R=Ce, Sm are contradictory concerning the na- 
ture of the phase transition from magnetic order to super- 
conductivity as a function of doping. The phase diagram 
of the Ce system reported by Zhao et al. [4] allows both, a 
magnetic and a possibly structural QCP, whereas Sanna 
et al. [7] and Shiroka et al. [8] claim "nanoscopic" coex- 
istence of magnetic order and superconductivity (phase 
separation with nanometer sized domains) which indi- 
cates the absence of a magnetic QCP (due to a first or- 
der phase transition as a function of doping). The lat- 
ter situation agrees with the claim of (nanoscopic) phase 
separation between magnetism and superconductivity in 
the Sm system[9]. Alternatively the /xSR data presented 
by Sanna, and Shiroka and co-workers [7-9] may be ex- 
plained by microscopic coexistence of magnetic order and 
superconductivity, i.e., they coexist (spatially) homoge- 
neous and form a single thermodynamic phase that is 
both superconducting and magnetically ordered. This 
implies that both order parameters are coupled, or, in 
different words, the same electrons contribute to super- 
conductivity and magnetic order[10]. 

The structural phase transition for i?=Sm is discussed 
controversially: Margadonna et al. [11] reported a struc- 
tural phase transition with a strongly suppressed phase 
transition temperature Ts in underdoped superconduct- 
ing SmFcAsOi_a;Fa;, whereas Martinelli et al. [12] find 
that Tg is almost doping independent even beyond the 
optimally doped regime. In the present work we want 
to clarify these contradicting results with new experi- 
mental data and a comparison with available literature. 
We will present a thorough analysis of our experimental 
data that reveal, new details of the phase diagrams of the 
i?FeAs0i_2:Fa; (i?=Sm, Ce) superconductors. 

The intricate interplay between structure, magnetism 
and superconductivity requires the combination of sev- 
eral experimental techniques. No single experiment ex- 
ists that can detect all relevant aspects. It is of utmost 
importance to realize the limitations and advantages of 
each experimental technique. For example, diffraction 
experiments with neutrons or X-rays are only sensitive 
to order with long correlation lengths. This is an es- 
sential limitation. For example, when discussing micro- 



scopic coexistence and phase separation of SC and small 
moment, or short range magnetic order. Local probes 
such as Mossbauer spectroscopy and muon spin relax- 
ation (/xSR) on the other hand, are very sensitive to dis- 
ordered magnetism. /xSR is also an excellent tool to de- 
termine the absolute values of the magnetic penetration 
depth of superconductors. 

Here we report the results of our synchrotron X- 
ray diffraction (XRD), muon spin relaxation and ro- 
tation (/^SR), Mossbauer spectroscopy, and electri- 
cal resistivity experiments of CcFeAsOi-ajFj; for x=0, 
0.042(2), 0.048(3), 0.063(2), 0.145(20), 0.150(20), and 
SmFeAsOi_^F^ with x=0, 0.02, 0.04, 0.06, 0.08, 0.10. 
We also present specific heat measurements of the un- 
doped Sm system and of a;=0.06. The results we ob- 
tain for the Ce system clarify the above mentioned fea- 
tures of the phase diagram. We unambiguously show 
that magnetic order and the tetragonal to orthorhombic 
structural transition are suppressed with doping but only 
short range magnetic order (SRO) is found together with 
superconductivity for a;=0. 063(2) in the Ce, and x=:0.06 
in the Sm system, whereas the structural transition dis- 
appears as superconductivity emerges with doping. 

Our experiments confirm a mixed phase of supercon- 
ductivity and magnetic order in both systems but the 
indications for phase separation are not sufficient to 
rule out a magnetic quantum critical point that possi- 
bly emerges as a function of F doping. We determined 
the structural phase transition temperature from a de- 
tailed analysis of the lattice distortion of the orthorhom- 
bic phase. In several publications [4, 12] the broadening 
of Bragg peaks was used to obtain the structural transi- 
tion temperature. Our results show that this procedure 
can be inaccurate in the systems at hand. In fact, this 
kind of analysis should always be supported by additional 
experimental evidence such as thermal expansion or spe- 
cific heat to verify the presence of a phase transition. 
To study the order parameter of superconductivity we 
employed transverse field /LtSR. We find that the temper- 
ature dependence of the superfluid density ti^cxI/A^ is 
consistent with a nodeless s-wave gap in both systems. 
This is in agreement with measurements of the penetra- 
tion depth A that find an exponential decrease as a func- 
tion of temperature [13, 14]. We discuss our results in 
terms of the effects of random stress, and disorder on the 
magnetic, and structural phase transition. We also study 
the doping dependence of the rare earth magnetic order. 

EXPERIMENTAL METHODS 

Polycrystalline samples of SmFeAsOi-xFj, and 
CeFeAs0i_a;F2: were synthesized using the two step 
solid state reaction route described by Kondrat et al. 
[15] and annealed all samples in vacuum. The crystal 
structure was verified by powder Mo Ka X-ray diffrac- 
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tion [11, 15-18] and the FullProf software[19]. Electrical 
resistivity was measured using a standard 4-point 
geometry with an alternating dc current. Synchrotron 
powder X-ray diffraction data was collected (A=l. 07813 
A) for the (2,2,0)t and (0,0, 6)t Bragg peaks (subscript 
T refers to tetragonal), in a temperature range between 
7 K and 200 K, at the beam line for Resonant Magnetic 
Scattering and High-Resolution Diffraction (MAGS) at 
the 7 T Wiggler at BESSY in Berlin, Germany. Muon 
spin relaxation and rotation (^SR) experiments were 
conducted at the Paul Scherrer Institute (PSI) using 
standard *He flow cryostats, a '^He-^He dilution cryostat 
at the GPS; Dolly, and LTF instruments. Mossbauer 
spectroscopy was conducted using a ^Hc flow cryostat 
with a '^^Co-in-Rh matrix gamma radiation source kept 
at room temperature (emission line half-width-at-half- 
maximum of 0.130(2) mm/s). From X-ray diffraction, 
less than 5 mol% impurity phases are inferred. From 
Mossbauer spectroscopy less than 1 atom% of Fe bearing 
impurities are inferred for CeFeAsOi-ajF^; with x=0, 
0.48(2), and 0.63(3); x=0.145(20) contains «5 weight% 
of Fc2As. Panarina et al. [20] studied samples from the 
same batch of SmFcAsOi-ajFj, for x>0.06 and reported 
evidence for bulk superconductivity from a carefully 
analysis of susceptibility measurements. We determined 
the doping levels of the Sm samples using wavelength 
dispersive X-ray spectroscopy. To determine the F 
contents of the Ce samples we used nuclear quadrupole 
resonance, with the corresponding study of the electronic 
properties to be published separately [21]. 



SmFeAsOi_:,;F^ 

In this section we present the results of our electrical 
resistivity (Sec. ), synchrotron X-ray diffraction (Sec. ), 
and muon spin relaxation measurements (Sec. ). We also 
present specific heat measurements in Sec. , and to sup- 
port the results of our /iSK measurements. The struc- 
tural transition and the Fe magnetic transition are both 
broadened. /iSR enables us to disentangle the temper- 
ature dependence of the magnetic order parameter and 
the magnetically ordered volume fraction — this is, in gen- 
eral, not possible with diffraction experiments. In Sec. 
we use this knowledge to construct a phcnomenological 
model that correlates all static experimental parameters 
of the /xSR in a self consistent way. In Sec. we discuss 
the evolution of the Sm magnetic order as a function of 
doping and in Sec. the magnetic fluctuations unique to 
SmFeAsOi-aiF-c. In Sec. we summarize our main re- 
sults and present our phase diagram of SmFeAsOi_^F^. 
Phase transition temperatures are shown in Tab. I, p. 19. 




Temperature (K) Temperature (K) 

FIG. 1. Normalized electrical resistivity as a function of tem- 
perature of SmFeAsOi-a;Fii. (Blue) Red curves refer to (non) 
superconducting compounds. Curves are shifted vertically for 
clarity. 



Electrical resistivity measurements 

We discussed the temperature dependence of the elec- 
trical resistivity p{T) of SmFeAsOi-a^F^: in Fig. 1 in 
detail in a previous work [22]. Here we briefly men- 
tion the essential features for later comparison with the 
CeFeAsOi-a^F^; system. For x=0, and 0.02 we observe a 
cusp in the resistivity that has been used in the past to 
determine the structural and magnetic phase transition 
temperatures at the maximum and the inflection point, 
respectively [22, 23]. We find the maximum of the resis- 
tivity at rP„^=160, 146, and 135 K for x=0, 0.02, and 
0.04, respectively (see Ref. 22 for details). For the inflec- 
tion point we flnd i;''=136 K, 125 K, 97 K for x=0, 0.02, 
and 0.04, respectively. Upon doping, this anomaly is 
slightly suppressed to lower temperatures. It is absent in 
the superconducting compounds (a;>0.06). Instead, p{T) 
is linear for r>150 K and exhibits a pronounced change 
of slope at 150 K (see Fig. l(b)-(d)) which has been 
attributed as to the spin fluctuations related to proxi- 
mate spin density wave state. We find Tc=36.2, 44.5, 
and 52.1 K for a;=0.06, 0.08, and 0.10, respectively. 



X-ray diffraction measurements 

To study the tetragonal to orthorhombic structural 
transition, we measured the tetragonal (2,2,0)t Bragg 
peak (space group PA/nmm) as a function of tem- 
perature. In Fig. 2 we show representative diffrac- 
tion patterns of SmFeAsO. When we lower the tem- 
perature below the structural transition temperature 
Ts=143.7(-Hl-5) K, 134.5(-Hl-5) K, 106.3(-Kl-5) K for 
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FIG. 2. Left: representative X-ray diflFracted patterns of the (2,2,0)t tetragonal and the (4,0,0)o and (0,4,0)o orthorhombic 
Bragg peaks at different temperatures for SmFeAsO. Right: orthorhombic distortion S{T) as a function of temperature. Colored 
symbols indicate the temperature range in which S{T) is concave, whereas the gray symbols indicate that 5{T) is convex. Inset: 
5(T — >■ 0) as a function of the bulk critical temperature Ts from this work and from Margadonna et al. [11] 
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FIG. 3. Lattice constants as a function of temperature for 
SmFeAsOi-xFa:. The black circles refer to the a and b lattice 
parameters and the red triangles to c. Below the tetragonal 
to orthorhombic transition, a and b constants are divided by 
V2 for comparison. 



x=0, 0.02, and 0.04, respectively, the (2,2,0)t peak splits 
into the (4,0,0)o and (0,4,0)o orthorhombic Bragg peaks 
(space group Cmma). 

The order parameter of this transition is the or- 
thorhombic distortion S{T) in Fig. 2. It is the normalized 
splitting of the Bragg peaks S{T)={a — b)/{a + b) where 
a and h are the lattice constants. We determined the 
temperature dependence of the lattice constants shown 
in Fig. 3 from the Bragg peak positions. The struc- 
tural transition of a::=0, 0.02, and 0.04 is indicated by 
a{T)^b{T) below the transition temperature Ts- We 
observe no such feature for a;=0.06, and 0.10 (we did 
not study x~0.08 with synchrotron X-ray diffraction). 
The temperature dependence of the specific heat Cp{T) 



of a;=0.06 in Fig. 11, p. 12 also shows no indication of a 
structural transition. This indicates that for a;>0.06 the 
structure remains tetragonal down to the lowest temper- 
ature T«10 K, that we measured in this study. 
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FIG. 4. FWHM of the (2,2,0)t Bragg peak (left) above the 
structural transition, and (right) in the superconducting su- 
perconducting regime of SmFeAsOi-iF^;. Arrows indicate 
the superconducting transition temperature Tc, and the on- 
set of magnetic short range order at Tsro- Please note the 
smaU increase of the FWHM for 2;=0.10 below 200 K that is 
followed by a plateau between 150 K and ~100 K (see text 
for a detailed discussion of this anomaly). 

In the following we present a detailed analysis of the 
structural transition. When we approach the structural 
transition from high temperatures the tetragonal (2,2,0)t 
Bragg peak broadens (see Fig. 4) by up to ~70% before 
we can observe a clear splitting at Ts shown in Fig. 2. A 
priori it is unknown whether this is due to an unresolved 
splitting or an intrinsic broadening of the (2,2,0)t Bragg 
peak or a combination of both. When we approach the 
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structural transition from low temperatures the splitting 
of the orthorhombic (4,0,0)o and (0,4,0)o Bragg peaks 
decreases and becomes comparable to the full-width-at- 
half-maximum (FWHM) of the Bragg peaks. In this 
temperature range the Bragg peak positions and their 
FWHM become correlated fit parameters. To lift this 
correlation we keep the FWHM at the constant value wq . 
The value of S{T) for T>Ts close to the transition de- 
pends on the choice of wq wheras at lower temperatures 
T<Ts, S{T) is independent of wq- Here we set wq to 
the value that we obtain for TwlO K. For the discussion 
of the order parameter we therefore limit the analysis to 
data for which the splitting is larger than wq. This coin- 
cides with the temperature range for which a description 
of the data with a single Bragg peak is not sufficient any- 
more. 

The so obtained orthorhombic distortion is shown in 
Fig. 2. In a first step we attempt to analyze it with a gen- 
eral power law 5{T)(x{l — {T /Ts)")^ ■ Here a controls the 
temperature dependence of the order parameter at low 
temperatures and the critical exponent (3 the steepness 
close to Tg. Below Tg the bulk orthorhombic distortion 
follows the above power law. However, this analysis docs 
not yield reliable estimates of the critical exponents. 
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FIG. 5. The normalized orthorhombic distortion of 
SmFeAsOi-a:Fa: compared with the temperature dependence 
of a 2D-Ising, and 3D-Ising type order parameter. Close to 
the transition it is not clear whether 5{T) is the true bulk 
order parameter. Therefore the temperature error is large 
and asymmetric (see text). Qualitatively all three order pa- 
rameters are closer to the 2D-Ising than the 3D-Ising type. 
However, 5{T) tends towards the 3D-Ising type but within 
the error we can make no clear classification. 

Nevertheless, we are able to identify trends of the order 
parameter graphically. We normalize the order param- 
eter to the value at T=0.6Ts and the temperature axis 
with the estimated transition temperature Ts in Fig. 5 
and compare it with the 2D-Ising S{T)(x{l — T/Ts)'^ with 
/3=0.125, and the 3D-Ising (3=0.325. In this analysis we 
have to take into account the large uncertainty of T5. For 



all three doping levels, the order parameter is close to 
the 2D-Ising universality class. However, 6{T) of x=0.04 
tends towards the 3D-Ising case. A clear classification is 
not possible due to the large error of T5. 

This result is similar to that obtained by Wilson et al. 
[24] for the magnetic transition. They concluded a 
change of dimensionality from 2D- to 3D-Ising of the 
magnetic phase transition under doping by an analysis 
of available experimental data. We find a similar trend 
for the structural transition of the Sm system. In sec- 
tion , p. 19, we show that this trend is more pronounced 
in the Ce system (due to the slightly higher doping level) . 

In the inset of Fig. 2 we show (5(7— >0) as a function 
of Ts from our, and the work by Margadonna et al. 
[11]. We find, that for low doping levels for which 
r5>100 K, (5(T-^0) depends linearly on Ts and a fit 
yields 5{Ts)={0.02A{l)K-^Ts - 0.8(1)) • IQ-^. The two 
data points with Ts<70 K in the inset of Fig. 2 that 
deviate from this linear dependence are from supercon- 
ducting materials studied by Margadonna et al. [11] This 
indicates that in a narrow part of the superconducting 
regime of SmFeAsOi-ajF^; the orthorhombic distortion is 
stabilized while the critical temperature is suppressed by 
doping. 
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FIG. 6. The relative change of the lattice constants in the 
proximity of Tc of SmFeAsOi-aiFj; with a;=0.10. The change 
of slope from negative {T<Tc) to positive {T>Tc) indicates 
that Tc decreases under pressure. The steeper slope of c{T) 
suggests that Tc is more sensitive to uniaxial pressure along 
the c-axis than to in-plane pressure. We were not able to 
resolve the effect of superconductivity for 2;=0.06 — close to 
Tc, its lattice parameters are constant within the experimental 
error (not shown). 



For .T=0.10 we observe superconductivity with 
Tc=52.1 K (from electrical resistivity). The structural 
transition is absent for this composition. But we find 
an upturn of the a, and c lattice constants below TfnTc. 
This upturn is shown in Fig. 6. The relative changes of 
the lattice constants at Tk,Tc with respect to the value at 
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r=10 K are Ac(Tc)/co<c(Tc)/c(10 K)-l=-1.7(2)•10-^ 
and Aa(Tc)/ao<a(Tc)/a(10 K)-l=-0.15(2)-10-'*. 

The origin of this anomaly of the temperature depen- 
dence of the lattice constants may be either magnetic or- 
der or superconductivity. But for a;=0.10 we can rule out 
a magnetic origin because our ^SR measurements show 
the absence of any magnetic order in this temperature 
range. It is more likely that the pressure dependence 
of Tc governs the lattice contraction. According to an 
Ehrenfest relation, the jump Aa^ of the thermal expan- 
sion coefficient ai at the derivative of the thermal 
expansion along a crystallographic direction i, is propor- 
tional to the pressure dependence dTc/dpi of Tc (see e.g. 
Ref. 25). The contraction of the lattice upon heating to- 
wards Tc and subsequent expansion suggests a jump from 
a<0 below Tc to a>0 above Tc, AaocACpdTc/dp should 
then be negative. Since c(T) is steeper than a(T), the 
jump AacOcdTc/dpc is larger than for the a-axis. This in- 
dicates that the pressure dependence of Tc is dominated 
by uniaxial pressure along the c-axis. This anisotropy of 
the pressure dependence should be kept in mind when 
comparing pressure experiments with varying degrees of 
hydrostaticity and different pressure cell geometries [26]. 

Qualitatively, this result agrees with the decrease of Tc 
of nearly optimally doped SmFeAsOi^ajF^; under high 
pressure[27]. Quantitatively, c(Tc)/co<-1.7(2)-10"'*, 
and Aa(Tc)/ao<— 0.15(2)-10^'* are consistent with 
the results of thermal expansion experiments of 
BaFe2_a;Coa:As2 by Bud'ko et al. [28]. 

For a;=0.06, we were not able to resolve any effect on 
the lattice at or below T^. The jump of the thermal ex- 
pansion coefficients AaccACpdTc/dp is proportional to 
the specific heat jump at the transition, which should 
be proportional to Tc [29]. This suggests a 30% smaller 
effect on the lattice constants of a:=0.06 compared to 
x=0.10 based solely on the different Tc- In addition, ab- 
solute value of dTc/ dp~-2 to -6 K/GPa for nearly opti- 
mally doped SmFcAs0i_2;Fa; and SmFeAsOi_5 [30-32] is 
reduced by a factor of approximately two to dTd dp~+l 
to +2.Q K/GPa for underdoped SmFcAsOi_2;F^[30, 33]. 
These two factors equate to a reduction of the effect on 
the lattice constants (the jump Aa) by w70%. Over- 
all, this should make the effect of superconductivity on 
the lattice constants of a;=0.06 undetectable in our X-ray 
diffraction study. 

Fig. 4 shows the fuU-width-at-half-maximum (FWHM) 
of the (2,2,0)t peak for a;=0.06. When cooling, it in- 
creases continuously until Tc and then decreases slightly. 
For 1=0.10, it is reduced by 50-70% compared to a;=0.06 
and increases continuously without any sharp anomalies. 

The FWHM is a measure of the structural disorder 
and its inverse is related to the correlation length of the 
structural long range order. For the sample with a;=0.06, 
lattice fluctuations are enhanced upon cooling until we 
reach Tc- In contrast, structural fluctuations are signif- 
icantly reduced for x—0.\Q but also show a similar in- 



crease towards lower temperatures, but no anomaly at Tc- 
The behavior of the FWHM in figure 4, is reminiscent of 
the evidence of spin fluctuations observed in supercon- 
ducting LaFeAsOi_2;Fa; by nuclear magnetic resonance 
[34, 35]. The importance of magneto-elastic coupling has 
been discussed also in other theoretical and experimen- 
tal studies [36-39]. This should play an important role 
in the temperature dependence of the FWHM. Tropeano 
et al. [40] and Zhu et al. [41] suggested that magnetic 
fluctuations observed above Tc are suppressed in favor of 
superconductivity. 

A suppression of spin fluctuations by superconductiv- 
ity may also cause a decrease of the FWHM below Tc- 
We observe a decrease of the FWHM below Tc only for 
a;=:0.06. However, for x=0.10, the system may be suffi- 
ciently far away from a spin density wave instability and 
for this mechanism to be relevant 

Alternatively, stresses and disorder may similarly 
broaden the structural transition and cause a temper- 
ature dependence of the width of the tetragonal Bragg 
peak in a wide temperature range [42-46]. 

Margadonna et al. [11] reported that superconduc- 
tivity can be found in the orthorhombic phase of 
SmFeAsOi_2;F^ in a narrow doping range. In Fig. 24 
we show the orthorhombic distortion as a function of 
the structural transition temperature from this work and 
from Margadonna et al. [11]. It is evident that both 
datasets are consistent. 

In contrast, Martinelli et al. [12] reported evidence 
for a tetragonal to orthorhombic structural transition 
with r5sal70 K in superconducting SmFeAsOi-ajF^, with 
x—Q.lQ . They obtained this result by a microstrain 
analysis that showed an anisotropic broadening of the 
(1,0, 0)t Bragg peak. They concluded an unresolved 
splitting of the Bragg peak due to the structural tran- 
sition. This study was limited to temperatures above 
90 K. The authors suggested that at temperatures lower 
than 90 K a splitting of the Bragg peak should appear. 
We also observed a significant broadening of the tetrago- 
nal (2,2,0)t Bragg peak of the superconducting materials 
over a wide temperature range, as shown in Fig. 4. The 
Bragg peak FWHM for a;=0.10 even shows a plateau like 
feature between 90 K and 150 K with a small drop of 
the FWHM at higher temperatures. This resembles the 
anomaly observed by Martinelli et al. [12] that they in- 
terpreted as an indication for a structural transition. But 
we find no splitting of the Bragg peak down to «10 K 
that would indicate a structural transition. Considering 
that F doping suppresses Ts to 107 K already for a;=:0.06 
it seems unlikely that for a;=0.10 the transition should 
occur at 170 K as suggested by Martinelli et al. [12]. 
Here we point out that in the Fe based superconductors 
an increase of the FWHM of a Bragg peak alone cannot 
be taken as sufficient evidence for a structural transi- 
tion. This becomes even more evident when comparing 
the structural transition temperature obtained from the 
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analysis of the FWHM for the non superconducting ma- 
terials. With this analysis Martinelli et al. [12] find only 
a small decrease of T5 from ^175 K in the undoped ma- 
terial to sal70 K for a;=0.1 whereas we clearly observe a 
significant suppression of Ts from «144 K in the undoped 
material to «105 K for x=0.06 by the analysis of the or- 
thorhombic distortion. Please take note that the analysis 
of the FWHM in Ref. 12 yields transition temperatures 
that are up to «60 K too high. 

Muon spin relaxation measurements 

Introduction 

For the fiSK experiments a continuous beam of S'=l/2, 
nearly 100% spin polarized positve muons is directed onto 
the sample. Through inelastic Coulomb scattering with 
the sample the muons completely lose their kinetic en- 
ergy and finally stop at interstitial lattice sites. They will 
usually stop in the sample at depths of 0.2-0.3 mm. At 
the final interstitial lattice site the muon spin processes 
in the local magnetic field. For a single measurement 
several million positrons, generated by the muon decay, 
are recorded in time histograms by a set of positron de- 
tectors. Properly normalized, these histograms yield the 
time evolution of the muon spin polarization P(t). This 
is possible because, during the muon decay the positron 
is emitted preferentially along the direction of the muon 
spin. The analysis of P{t) in zero, longitudinal, or trans- 
verse magnetic fields can reveal properties of, e.g. mag- 
netic order, magnetic fluctuations, and superconductiv- 
ity, respectively. 

In a previous work we calculated possible muon sites. 
In agreement with experiment we found two muon sites 
(see below and Ref. 47). In a recent DFT study, De Renzi 
et al. [48] calculated similar muon sites. Accordingly we 
chose a model for commensurate magnetic order to de- 
scribe the zero field ^SR time spectra that we deduced 
from the following general model with the two muon sites 
A and B: 

P{t) -K„a<,(T)[a(^e-^^-^*cos27r/^t + ie-^^-^*) 

+ (1 - a)(^e-^--«* cos 27r/si + ^er^^-^*)] 

+ (1 - V,nagiT))[a{Git,aA)e-^""^-^ 
+ il-a)G{t,aB)e-^—'^]. 

The first term models the signal from muons that come 
to rest in the volume of the sample that is magneti- 
cally ordered. Its volume fraction is the magnetic vol- 
ume fraction VmagiT). The second term is due to muons 
that stop in the volume of the sample that does not 
show long range magnetic order. Its volume fraction 
is {l—Vmag{T)). Each term contains two muon sites A 



and B with average muon occupation a and (1 — a), 
respectively. The last term models the paramagnetic 
part of the sample were muon-spin relaxation is caused 
by nuclear magnetic dipoles, described by the Gaussian- 
Kubo-Toyabe relaxation function G(t,a), and electron 
spin fluctuations, described by the exponential term [49]. 

The model (1) contains many free parameters - not all 
can be determined accurately. It is often impossible to 
distinguish between the dynamic relaxation rates Al and 
Xnm in the magnetically ordered, and the paramagnetic 
phase. Therefore, we set XL,A=Xnm,A and XL,B=Xnm,B- 
The same applies to the rates that characterize the nu- 
clear dipole field distribution, i.e., we set a=aA=o'B- In 
case of a;=0.06 we only observe one strongly damped sig- 
nal without a frequency but two dynamic rates A^.^i and 
Xl.b we cannot distinguish At, a, and At,s, therefore we 
set Xt=Xt,a=Xt,b- For x=0, 0.02, and 0.04, above the 
rare earth magnetic ordering temperature for the Sm sys- 
tem we never observed more than one frequency and we 
set a=l. In summary model (1) reduces to 

P{t) = V„,agiT)[^e~^^'cOs2TTf^t 

3 

+ (1 - V,nag{T))G{t, a) [ae-^^'^ + {I - a)e-^^'«] . 

(2) 

Above model describes the signal of the majority phase 
in the sample. Small amounts of magnetic impurity 
phases cause relaxation of a few percent of the total sig- 
nal during the dead time of the detectors for a;>0.06, this 
reduces P{t = 0) below 1. 

In Ref. 47 we have presented a calculation of the muon 
spin precession frequency based on our muon site cal- 
culation and, at the time available estimates of the Fe 
ordered moment from neutron diffraction. For LaFeAsO, 
de la Cruz et al. [18] reported an Fe moment of psO.36 /jb- 
Using this value we calculated the dipole fields Bioc and 
the muon frequencies /=7^/(27r)_B/oc at the two muon 
sites. However, this magnetic moment yields frequen- 
cies that are by a factor of 1.86 smaller than the ex- 
perimentally reported frequencies of 7^=23 MHz and 
/b=3 MHz [6, 23]. Only recently, Qureshi et al. [50] 
reported a neutron diffraction study of LaFcAsO with 
different samples. They find a much larger Fe mag- 
netic moment of 0.63(l)/iB. Using this value, the muon 
frequencies we calculate increase by factor of 1.75 to 
/a=21.7(3) MHz and /s=2.3(4) MHz where the error is 
estimated from the error of the magnetic moment given 
by Qureshi et al. [50]. Within this error, the calculated 
and measured muon frequencies deviate by only wlO- 
20%. 

This result has several implications: 

1. The coordinates of the calculated muon sites are 
very close to the real muon sites. 
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FIG. 7. Panels (a) through (f) show the zero field /.iSR time spectra characteristic for the SmFeAsOi-^Fa; compounds we study 
here. For a;<0.04 we observe coherent muon spin precession which proofs long range magnetic order. For a;=0.06 a growing 
fraction of the signal shows fast relaxation below Tsr.o=28(2) K which indicates short range magnetic order. Panel (g) shows 
longitudinal field /iSR measurements for a;=0.10. Small magnetic fields of up to 100 mT do not have any effect on the time 
spectra. This indicates that only time dependent magnetic fields, i.e., spin fluctuations cause the relaxation and not static 
magnetic moments. 



2. The interaction of the muon with its surrounding, 
i.e., electrostatic polarization, magnetic polariza- 
tion, and local deformation of the crystal lattice 
have negligible effects on the /zSR measurements of 
the studied ferropnictides, i.e., with ^SR the intrin- 
sic magnetic properties are measured. 

3. No Fermi-contact hyperfine field contributes to the 
local magnetic fields at the interstitial muon lattice 
sites, because we can explain them fully by mag- 
netic dipole fields. 

These conclusions are in broad agreement with the re- 
sults of De Renzi et al. [48]: Based on their calculation 
of the muon site, they determined a Fe magnetic moment 
of 0.68 fiB, which is only 10% larger than the value deter- 
mined by neutron scattering [50]. They also conclude that 
Fermi contact hyperfine fields, by symmetry, do not con- 
tribute to the local magnetic fields at the muon sites [48]. 
In addition they find a third muon site that may be rel- 
evant for muon diffusion, which takes place at tempera- 
tures larger than, for ferropnictides typical magnetic and 
structural phase transition temperatures[51]. 

The magnetically ordered phase may have a magnetic 
correlation length that is shorter than a few nanometers 
or may be broken up into small and isolated clusters. In 



theses cases the relaxation rate At can be much larger 
than the precession frequency As a result only a fast 
decay of the muon polarization without coherent muon 
spin precession can be observed. 

For SmFeAsOi-i^F^; we observed both situations, long 
range order with, and short range order without muon 
spin precession: a well defined muon spin precession, 
shown in Fig. 7 due to long range magnetic order for 
x<0.04, and for a;=0.06 only a fast relaxation indicat- 
ing short range magnetic order. For higher doping levels 
static magnetic order is absent. In the following we will 
discuss the results of the measurements in detail. 



Fe magnetic order 

In Fig. 7 we show our /iSR time histograms of the Sm 
compounds. From the data we deduce, that the magnetic 
transition in the non-superconducting SmFeAsOi_a;Fa; 
with x=0, 0.02, 0.04 occurs in two steps. We illustrate 
this process for SmFeAsO in Fig. 8 with a detailed series 
of histograms at different temperatures: 

Firstly, the data display coherent muon spin preces- 
sion already at high temperatures where only 5-10% of 
the volume are magnetically ordered. In the temperature 
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range between wlSO K and «140 K the magnetic volume 
fraction Vmag{T) in Fig. 9 increases gradually. However, 
the temperature dependence of the order parameter, the 
muon spin precession frequency /^(T), shown in Fig. 10 
does not follow the typical temperature dependence of 
an order parameter /^(T)oc(l — {T/Tm)"'Y — this is indi- 
cated by the symbols colored in gray in Fig. 10. Instead, 
ffj,{T) remains almost constant or decreases with decreas- 
ing temperature. As the precession amplitude increases, 
the damping rate Xt(T) of the muon spin precession in- 
creases to w20 MHz, see Fig. 10 and Fig. 8. 

Secondly, only when 80-90% of the sample volume are 
ordered, the damping rate Xt{T) decreases and ffi{T) 
displays the typical temperature dependence expected for 
an order parameter /^(T)cx(l - (T/Tn)")'^- 




time t (\is) 



FIG. 8. Zero field muon spin relaxation time spectra of Sm- 
FeAsO (curves are shifted vertically for clarity). We find 
muon spin precession, the indicator for long range magnetic 
order, already at temperatures for which only 5-10% of the 
material show magnetic long range order. This indicates that 
more and more clusters become ordered as we lower the tem- 
perature. Only close to percolation f^{T) follows the typical 
temperature dependence of an order parameter (see Fig. 10). 

Qualitatively the transition follows the same pattern 
for a:=0, 0.02, and 0.04: A broad magnetic transition 
with a gradual increase of Vmag{T) (see Fig. 9) with an 
unusual temperature dependence of fp,{T) and a peak in 
\t{T) (see Fig. 10). We will present a phenomenological 
interpretation of this broad transition in the next section 
(see p. 11). 

Wc also find a broad magnetic transition for a;=0.06. 
magnetic fluctuations of the Sm 4/ moments (see Sec. , 
p. 13). At high temperatures, the relaxation of P{t) is 
due to the magnetic fluctuations of the Sm 4/ moments 
(see Sec. , p. 13). The magnetic order occurs below 
Tsro=33(l) and causes, as we lower the temperature, a 
gradually increasing part of P{t) to relax during the first 
0.2 /is of the time histograms in Fig. 7. This is indicated 
by the gradual increase of Vmag{T) that reaches «100% 



below Tioo=6(2) K, and At(T) below T^ro (see Figs. 9, 
and 10). 

However, down to a temperature of 1.8 K we find 
no coherent muon spin precession — an indication for 
short range magnetic order. Since the susceptibility 
measurements [20] also show bulk superconductivity, we 
conclude that below Tioo the material is 100% magneti- 
cally ordered and 100% superconducting. Sanna, Shiroka 
and co-workers [7-9] take this as evidence for nanoscopic 
coexistence, i.e., phase separation on length scales of a 
few nanometers. However, this would only be clear ev- 
idence for phase separation if the /uSR signals (volume 
fractions) from both phases were clearly separable, as 
e.g. in Bai_2;K2;Fe2As2 [52-55]. Clear evidence for micro- 
scopic coexistences would be coupled order parameters as 
e.g. in BaFe2-a;Coa;As2[56]. Considering the available ev- 
idence for the present case, we suggest that the situation 
is consistent with both situations. 

At higher doping levels a;=0.08, and 0.10 we observe 
no sign of Fe magnetic order. Instead, the relaxation of 
P{t) in zero field is mostly due to magnetic fluctuations 
of the Sm 4/ moments (see Sec. , p. 13). 

Considering this broad magnetic phase transitions, it 
is not straight forward to define a phase transition tem- 
perature. The temperature for which Vmag{T) becomes 
non-zero is T,™=151(1) K, 143(2) K, 115(5) K, 33(1) K 
for x~Q, 0.02, 0.04, and 0.06, respectively. We define 
Tioo as the temperature below which the majority of 
the sample is magnetically ordered and Vmag{T) sat- 
urates. We find rioo=130(2) K, 122(2) K, 83(1) K, 
6(2) K for a;=0, 0.02, 0.04, and 0.06, respectively. We 
choose the Neel temperature T/v so that f^ {T) follows a 
power law /^(r)oc(l- (T/Tjv)")^ for T<rjv=135.3(2) K, 
128.6(10) K, 86.7(22) K for a;=0, 0.02, and 0.04. 

The fits with this general power law are shown in 
Fig. 10. The exponent a=1.49(4), 3.0(3), 5.5(21) for 
a;=:0; 0.02, and 0.04, respectively, controls the satura- 
tion at low temperatures. The increase of a with doping 
indicates that the ordered magnetic moment saturates 
faster as a function of temperature the higher the F con- 
tent. The critical exponent /3=0.143(2), 0.19(2), 0.15(8) 
for a;=0, 0.02, and 0.04, respectively, controls the steep- 
ness of ffiiT) close to T/v and indicates the universality 
class of the phase transition. Due to the unusual tem- 
perature dependence of f^{T) in the vicinity of the 
uncertainty of /3 is too large to draw definite conclusions. 
The values of (3 are between /3=0.125 expected for a 2D- 
Ising transition and /3=0.325 of a 3D-Ising transition. 
This is consistent with previously reported values (see 
Ref. 24 and references therein). 

The saturation value of the precession frequency 
/o=23.07(4) MHz, 20.06(4) MHz, 11.91(7) MHz for x=0, 
0.02, and 0.04, respectively, is proportional to the ordered 
magnetic moment. In Ref. 47 we showed that the mag- 
netic moment of the i?FeAsO (i?=rare earth) is practi- 
cally independent of the rare earth. The ordered Fe mag- 
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temperature T (K) temperature T (K) temperature T (K) 

FIG. 9. The magnetic volume fraction Vmag with fits according to power law, and normal distributed transition temperatures 
(see Eq. (5)). A small but non-zero Vmag indicates magnetic order in small volumes of the sample that appear below Tsro- 
Vmag{T) Saturates at Tioo- For a-=0 we also show the result of a fit to simulated data (red triangles, see Sec. ). 




FIG. 10. Left: The temperature dependence of the muon spin precession frequency ,ffi,{T) of SmFeAsOi_a:Fi. It is proportional 
to the order parameter of the magnetic phase and indicates long range magnetic order. The gray symbols indicate that the 
precession signal is mainly due to isolated clusters and does not refiect the bulk order parameter. The black lines are fits to the 
data with the general power law described in the text — it was used to determine Tjv. The red arrows indicate the temperature 
of the maximum of At. Inset: The low temperature limit of the muon spin precession frequency ,fo=,ffi{T^O) as a function the 
Neel temperature Tjv. The line is a linear fit to the data /o(r;v)=0.224 MHz/K-riv-7.52 MHz. Right: The damping rate At 
of the muon spin precession signal is due to spatial and time dependent fiuctuations of the ordered magnetic moment as well 
as quenched disorder. However, within our model we can explain the peak at the transition with a large distribution of order 
parameters that narrows rapidly. The red lines shown in the left and right panel were generated with the phenomenological 
model that we present in the text. 



nctic moment of 0.61 of LaFcAsO[50] corresponds to 
a frequency of /o=23 MHz[6, 23]. This indicates that the 
measured frequencies in the Sm compounds correspond 
to an ordered Fe magnetic moment of 0.61 /is, 0.53 /is, 
and 0.32 for a;=0, 0.02, and 0.04, respectively. 

In the inset of Fig. 10 we plot /o = /^(T = 0) as a 
function of TV . We find that /o is linear in Tn and a fit 
to the data yields /o(rAr)=0.224 MHz/K-rjv-7.52 MHz. 
This seems to be a general feature of the ferropnictides: 
Fernandes et al. [57] analyzed magnetic neutron scat- 
tering data of BaFe2-a;Coa;As2 with a simple two band 
model and showed that, in agreement with experiment. 



the magnetic moment decreases nearly linear as a func- 
tion of the ordering temperature. However, it is not a 
general consequence of their two band model, but rather 
of the particular detuning of the Fermi surface pockets 
by doping[57]. 

The gradual increase of Vmag{T) suggests that mag- 
netic order occurs in clusters that show already the hall- 
mark of long range magnetic order: a non-zero muon 
spin precession frequency f^{T) shown in Fig. 10. The 
fact that close to Tat, f^{T) does not follow the typical 
temperature dependence of an order parameter indicates 
that the amount of ordered clusters has the main effect on 
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ffj,{T) and not the temperature dependence of the bulk 
order parameter. This is consistent with the peak in 
Xt{T). It suggests that near the transition the ordered 
clusters have a broad distribution of ordered moments, 
which dominates XriT) and narrows as the ordered mo- 
ments saturate at lower temperatures. In the next section 
we present a phenomenological model that quantifies the 
above considerations. 

This type of cluster formation could indicate a first 
order transition. However, specific heat measurements of 
the i?FeAsO type fcrropnictidcs show no clear sign of a 
delta-function like anomaly or a temperature hysteresis 
at the phase transition that would indicate a first order 
phase transition [58-60]. 



A phenomenological model for the Fe magnetic phase 
transition 



In the following we present a phenomenological model 
that will help to understand the broad transition that 
manifests in the slow increase of VmagiT), the peak of 
XxiT) and the unusual temperature dependence of f^(T) 
for T>Tn. 

The basis for this model is the following picture based 
on disorder. The investigated samples are polycrystals 
and each grain of a polycrystal may have a local Neel tem- 
perature Tj^ (as opposed to our estimation of the bulk 
T/v presented in the preceding section) which depends 
on density of defects and dopant atoms, and possibly the 
size and direction of stress. Therefore, we expect a prob- 
ability distribution P{Tj^) of local Neel temperatures: at 
a temperature T all grains with T'j^>T constitute the 
magnetic volume fraction: 



P{T'^)dT'^. 



(3) 



A priori it is not clear what functional form PiT'j^) 
should have. Authors of previous works chose the normal 
distribution[9, 24] 



P{T'n) = 



1 



(4) 



We show fits with a normal distribution and Eq. (3) in 
Fig. 9 as dashed lines. These fits systematically underes- 
timate the data close to T^ro- To improve the fit close to 
Tsro we construct a distribution based on a power law: 



6+1 



-(To - T)\ for To<T<To~ a2^^^''+^^ 



(5) 



0, else 



Here &>0 is an exponent and Tq the upper limit of 
the distribution. The prefactor is chosen so that a 



is a measure of the width of the distribution with 
Vmag{T=To—a)=0.5. The distribution is cut-off towards 
lower temperatures so that Vmag{T)<l. This distribu- 
tion yields the fits shown as gray lines in Fig. 9. It de- 
scribes the data well except close to Tioo, the temperature 
at which V„iag{T) saturates. The parameters we obtain 
from the fits are ro=167(2) K, 161(3) K, 139(4) K, and 
cr=29(2) K, 31(3) K, 44(4) K for x=0,0.02, and 0.04, 
respectively. Within the error 6=4.0(5) is independent 
of doping and was optimized simultaneously for all three 
doping levels. The exponent 6=4.0(5) may depend on the 
microscopic details of the disorder. It is close to 6=3.1 
predicted for the mean field McCoy- Wu model[61]. It de- 
scribes a disordered layered planar Ising magnet where 
the disorder is correlated perpendicular to the planes. 

In the following we will use the above power law distri- 
bution to also describe /^(T) and Xt{T). We do this by 
a simulation of zero field muon spin relaxation time spec- 
tra. P{t) is proportional to the real part of the Fourier 
transform of the probability distribution P{B) of the lo- 
cal magnetic field B at the muon stopping site: 



Pit) 



cos(7^B<)P(B)dB 



cosi2TTf^iT,Tj^)t)PiTjy)dT^. 



(6) 



Here, /^(r)=/o(r^)(l - (T/Tj^)")'' where a and /3 are 
variable fit parameters, and /o(T;;^r) =0.224 MRz/K-T^- 
7.52 MHz as discussed above and shown in Fig. 10. 

We fit the spectra simulated by Eq. (6) with a damped 
oscillation P{t, T)=Pq{T) cos(27r/^(r)t) exp(-AT(T)t) 
to obtain the frequency /^(T), the damping rate XriT), 
and the magnetic volume fraction which is just PoiT). 
We used this procedure to optimize the critical expo- 
nents a and (3 for temperatures lower than T/v=135.3 K. 
For .T=0 we find that a=1.35(3) and /3=0. 170(3) yield 
the best description of /^(T) for T<T7v=135.3 K. 
This fit is shown as a red line in Fig. 10. Although 
we only optimized the critical exponents for T<Tn 
the model gives good quantitative agreement also for 
T7v=135.3(2) K<r<T,ro=151(l) K. 

This procedure also yields quantitative agreement with 
relaxation rate At(T) shown in Fig. 10. In particular, it 
reproduces the peak of XxiT) near the phase transition. 
Note that a and /3 were only optimized with respect to 
/^(T). The quantitative agreement with At(T) indicates 
that this simple model captures most features of the mag- 
netic phase transition correctly. In particular it describes 
the unusual temperature dependence of f^{T) for T>T]s[ 
and the decay of Ar(T) for TkTn- 

In summary, we can describe f^i{T) and Xt{T) by a 
broad and asymmetric probability distribution of mag- 
netic ordering temperatures. We suggest that the forma- 
tion of an inhomogeneous magnetic state with paramag- 
netic and magnetically ordered domains is an intrinsic 
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temperature T (K) temperature T (K) 

FIG. 11. Left: The temperature dependence of the specific heat Cp{T)/T of SmFeAsO0.94F0.06 (upper panel). This material 
is superconducting and does not show the anomaly associated with the structural or magnetic phase transition [62]. Cp{T)/T 
shows a broad maximum at r«75 K similarly to LaFeAsOo.gFo.i. To extract the contribution of the Sm CEF levels we follow 
the procedure described by Baker et al. [62] and calculate the difference Acp{T)/T of the specific heat of SmFeAsOo.geFo.oe and 
LaFeAsOo.gFo.i (lower panel). Acp{T)/T displays a broad Schottky-type anomaly due to the thermal population of Sm CEF 
levels and a difference in Sommerfeld coefficients A7=20.3(4) mJ/(molK^). Below «50 K the onset of the superconductivity 
in the La system (T/"") and Sm system {T^"^), and the magnetic order in the Sm system limit this analysis. Right: Low 
temperature Cp{T) of SmFeAsO and SmFeAsOo.geFo.oe. The A-anomaly indicates a second order phase transition which is 
attributed to the Sm magnetic order. This anomaly is only slightly suppressed by doping. A small Schottky-type hump in the 
specific heat of SmFeAsO centered at 2 K (marked with the small, gray arrow) was also observed in Ref. 62 and is most likely 
related to an impurity. 



response of the magnetic system triggered by extrinsic 
disorder. Its systematic study by experimental and the- 
oretical methods may reveal important properties of the 
system without disorder. 
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FIG. 12. The temperature dependence of the damping rate At 
of the muon spin precession. Its increase at low temperatures 
marks the onset of the Sm magnetic order. Lines are guides 
to the eye. 



Sm magnetic order 

Additional muon spin precession signals in the un- 
doped material mark the onset of the Sm magnetic order 



below T^"""=4.66(l) K[47]. In the ordered phase of the 
Sm moments for x=0.02 we also observe additional spin 
precession signals, but the onset temperature is not as 
clear. 

The specific heat Cp{T) probes the onset of Sm mag- 
netic order very sensitively. We have measured Cp(T) for 
and 0.06 at low temperatures, see the right-hand 
side panel of Fig. 11. We find a A-anomaly with a max- 
imum at the transition temperature T^™''^^=5.2(1) K, 
and 4.3(1) K for a:=0, and 0.06, respectively. From a gen- 
eral point of view we can expected that T^'".A'<r^™.<=P 
because, although the Sm magnetic order parameter is 
non-zero for T<T^'^ ''P it may be undetectable with ^SR 
because (a) the order parameter is too small, (b) the local 
magnetic field is a vector sum of the fields created by the 
Fe, and Sm magnetic order, which reduces the effect of 
the Sm order on the muon spin precession frequency, (c) 
the (spatial) fluctuations of the order parameter destroy 
a coherent muon spin precession, or (d) the magnetic 
field of the ordered magnetic moment cancels out at the 
muon site. 

The A-anomaly is not accompanied by a divergence of 
the dynamic relaxation rate Al (T) , see Fig. 13. This indi- 
cates that the hyperfine coupling of the muon is too small 
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or the critical fluctuations are too fast to be detected 
by /kSR. Thus, T^"^ cannot be determined directly from 
Al(T) and its temperature dependence is not related to 
the critical fluctuations expected for this second order 
magnetic transition (a detailed discussion follows in the 
next section) . A posteriori we can identify T^™ with the 
increase of the damping rate Xt{T). In Fig. 12 we show 
Xt(T) for low temperatures. The increase of Xt{T) cor- 
rectly indicates the onset of the broadening of the spin 
precession frequency spectrum due to the Sm magnetic 
order. It coincides with the maximum of the A-anomaly 
of Cp{T). It enables us to identify the ordering temper- 
ature of the Sm moments using fiSK by the increase of 
Xt{T). We call the ordering temperature of the Sm mo- 
ments that we determine in this way T^™' =5.6(3) K, 
4.7(1) K, 4.7(1) K, and 4.4(2) K for 2^=0, 0.02, 0.04, and 
0.06, respectively. 

This approach is not possible for 2;=0.08, and 0.10 be- 
cause the Sm magnetic order causes no additional static 
relaxation. Instead the relaxation is fully due to spin 
fluctuations. We verified the dynamic character of the 
relaxation by measuring in a longitudinal magnetic field, 
see Fig. 7(g). The absence of any sizable decoupling indi- 
cates the dynamic nature of the relaxation. The increase 
of Xj^ f^A/B) below ftiS K could be caused by the Sm order 
(see next section and Fig. 13). Because this feature is 
present for all doping levels it is likely that Sm magnetic 
order persists also for a;=0.08, and 0.10 with an estimated 
T^™sa3.5(10)K. However, this only indicates a possible 
phase transition and should not be taken as proof for Sm 
magnetic order for a;=0.08, and 0.10. 

Spin fluctuations 

In the preceding sections wc discussed the Fc and 
the Sm magnetic order. In this section wc focus on 
the magnetic fluctuations characterized by a, Xl,a and 
Xl,b- X^ i^ji/B) are the muon spin-lattice relaxation rates. 
In SmFcAsOi_a;Fa: and CeFeAsOi_a;Fa; fluctuations the 
rare earth and the Fe magnetic moment could contribute 
to the spin-lattice relaxation rate. 

Several nuclear magnetic resonance (NMR) studies 
show an enhancement of magnetic fluctuations in the 
paramagnetic state above the Fe magnetic transition, and 
the superconducting transition [34, 35, 63, 64]. However, 
we, and other groups, find no such fluctuations in the 
i?FeAsOi_:,F^ materials with ^SR[6-8, 23, 51, 65, 66], 
indicating that the associated muon spin-lattice relax- 
ation rate is too small to be measured in the time window 
of ^SR (several fis) . 

For the i?FeAsOi_a;F3, materials with a magnetic rare 
earth i?=Ce, Pr, and Nd no sizable spin-lattice relax- 
ation in /xSR has been reported [7, 47, 66] — with the ex- 
ception of i?=Sm[5, 9, 47, 67, 68]. It is not immediately 
clear where these spin fluctuations originate. In super- 



conducting SmFeAsOi-ajFa; the Sm spin fluctuations lead 
to l/TiocT^^-S (from NMR[69]). The /iSR spin-lattice re- 
laxation rates Al.(a/s) contain additional contributions 
that are temperature activated, [9, 47, 67, 68] see Fig. 13. 
Drew et al. [68] suggested that they are related to spin 
fluctuations in the vicinity of the superconducting tran- 
sition. 

In Fig. 13, all materials show a temperature activated 
(step-like) temperature dependence of Xl^(a/b) close to 
rw50 K followed by a plateau or slight increase and a 
steep increase below «5 K. In materials with Fe magnetic 
order (x<0.06) an additional cusp is present in the vicin- 
ity of Tjq which is due to the enhanced spin fluctuations 
at the phase transition. Apart from this cusp, the tem- 
perature dependence of the spin-lattice relaxation rate is 
qualitatively independent of doping. This suggests that 
the enhancement of Xl^(^a/b) at «50 K is a general feature 
of all SmFeAsOi-xFa; compounds and related neither to 
superconductivity nor to Fe magnetic order but rather to 
Sm spin fluctuations. 

For a;=0, 0.02, and 0.04 the onset of the Sm mag- 
netic order is accompanied by the appearance of addi- 
tional precession frequencies, and two distinct longitu- 
dinal relaxation rates Xl,a and Al,b {Xl,a ~ Xl.b for 
T>T§"'). 

The second relaxation rate Xl,b for a;>0.06 for T > 
T^™ is most likely related to the population of the second 
muon site in the Sm-0 layer [47]. The muons that stop 
in the Sm-0 layer have a much larger hyperfine coupling 
to the Sm spin fluctuations which leads to an enhanced 
spin-lattice relaxation rate Xl,b<Xl,a- The probability 
for a muon to stop in the Fe-As layer is a in (2), and 
the probability to stop in the Sm-0 layer is 1-a. We 
find that a is temperature independent for T<150 K and 
a=0.37(l), 0.38(1), 0.46(1) for x=0.06, 0.08, and 0.10, 
respectively. 

A a double-logarithmic plot of Xl,b{T) for 
SmFeAsO0.92F0.08 is shown in Fig. 14. The tem- 
perature activated contribution (green dashed line) 
used previously [67] increases down to rw30 K and 
then saturates. Xl^b{T), however, continues to in- 
crease. For 6-20 K we can describe it with a power 
law AL,B(r)=8.9(4)T-0-42(2)_ r^^a^ exponent is con- 
sistent with a recent ^^F NMR work [69] . Therein, 
Prando et al. [69] report a ^^F spin-lattice relaxation 
rate l/TicxT-°-6(i). They interpreted l/Ti to be due 
to critical fluctuations close to the magnetic phase 
transition of Sm, which should lead to a maximum 
of 1/Ti for T = T^"". Clearly, Az„(a/b)(T) does not 
show a maximum at T = T^™ ^ 5 K (see Fig. 13). It 
follows that /iSR is not sensitive to the itinerant Sm 
spin fluctuations sensed by NMR, but rather to the local 
fluctuations, i.e., we have to consider the fluctuations of 
the Sm moment between different crystal electric held 
levels. 

According to Orbach[70], the temperature depen- 
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FIG. 13. The longitudinal relaxation rates of 

SmFeAsOi_a:Fa;. The open symbols mark the additional 
relaxation rates that we find for a;>0.6. The Fe magnetic 
order causes a small cusp at high temperatures. It is followed 
by a temperature activated contribution that saturates below 
~30 K and a power law increase down to Ri5 K. At even lower 
temperatures the Sm magnetic order could be responsible for 
the steep increase below Ri5 K. 



dence of the fluctuation rate between CEF levels (at 
high enough temperatures) due to two-phonon scatter- 
ing leads to a finite lifetime of the CEF ground state 
given by 



^phonon ^ jyph' 



onon ^ ^A/iksT) 



1, 



(7) 



where rP'^°"°" is the line width of the quasi-elastic CEF 
excitation, fc^ is the Boltzman constant. The transi- 
tion between the degenerate spin states proceeds via 
the first excited CEF level with the energy A (the so- 
called resonant- Raman or Orbach process). Orbach's 
treatment [70] neglects any kind of interaction of the rare 
earth moments with the conduction band electrons. In 
first order, interaction with the conduction band elec- 
trons causes a Korringa-type broadening of the quasi- 
elastic CEF excitations [71] T^'^-'^f cx l/r^'^-'^f cx T. 
Gotze and Schlottmann [71] considered corrections to 
the linear temperature dependence due to magnetic ex- 
change interaction (valid at low temperatures) that lead 
to F'^''"^-'' oc T", where a < 1 depends on the density 
of states at the Fermi level and the hybridization (valid 
for T 3> Tk, where Tk is the Kondo temperature). Such 
power laws are often found by NMR measurements of 
the nuclear spin-lattice relaxation rate in heavy-fermion 
compounds, e.g., in CeFePO [72]. 

Assuming that both processes are independent, the 
overall line width of the quasi-elastic CEF excitation is 
just the sum F = rP^°"°" + T^'^-^f , i.e.. 



Xl{T) oc t oc 



1 



^phonon ^Sd— 4/ 

1 

.A/{kBT) _ I + 



6T" 



(8) 



where b scales the two contributions to the quasi-elastic 
line width, and Al is the muon spin-lattice relaxation 
rate (assuming the fast fluctuation limit). 

Within the error, b ^ 10.8(4) • IQ--* is independent of 
doping, A and the exponent a is shown in Fig. 15. The 
two relaxation rates for a;>0.06 have (within error) iden- 
tical temperature dependencies, therefore we optimized 
a and A simultaneously for both rates. The data is well 
described by the fits of Eq. (8) (see Fig. 13), however on 
cooling below T « 5 K the data increases steeper than 
expected, and reaches a cusp for T « 3 K. As discussed 
in the previous section, this is not related to critical fluc- 
tuations close to the Sm magnetic transition. This devi- 
ation suggests that magnetic Sm-Sm interaction should 
be taken into account for a description of the data for 
T < 5 K (see e.g. Rcf. 73) — such an analysis is beyond 
the scope of this work. 

A determined above is consistent with the CEF 
levels determined from the specific heat Cp{T) (see 
Fig. 11): To determine the CEF excitation energies 
of SmFeAsO0.94F0.06 we used Cp(T) of LaFeAsOo.ogFo.i 
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FIG. 14. The temperature dependence of spin-lattice relax- 
ation rate \l,b{T) for SmFeAsO0.92F0.08- This rate belongs 
to the muon site in the Sm-O layer. The temperature ac- 
tivated behavior is indicated by the green dashed-dot line. 
Below ~3Q K as a power law, below ~5 K it increases steeper 
but not cis steep as expected for critical fluctuations at the 
Sm magnetic phase transition. 



lattice relaxations rates of muons that stop in the Fe-As 
layer {Xl,a) and the Sm-O layer {Xl,b)- Qualitatively, 
Xl.{A/b){T) are doping independent and mainly probe 
magnetic fluctuations of the Sm 4/ states but with dif- 
ferent hyperfine coupling strengths. This is supported by 
the quantitative agreement of the lowest CEF excitation 
Ai determined by /iSR and specific heat measurements 
in this, and Baker's study [62], see Fig. 15. Furthermore, 
the ^^F NMR studies by Prando et al. [69] show a power 
law divergence of l/Ti, that they attribute to the Sm 
4/ states, similar to the power law increase we find be- 
tween 5 K and 30 K, sec Fig. 14. However, we suggest 
that /iSR is sensitive to the local Sm fluctuations (due 
to phonons and magnetic exchange with the conduction 
band electrons), whereas "'^^F NMR senses the collective 
Sm spin fluctuations (due to Sm magnetic order at low 
temperatures) . The Sm order is not accompanied by crit- 
ical fluctuations in the time window of //SR. The upturn 
of ^L.(A/B)(^^ < 5 K) should be due to the broaden- 
ing/splitting of the CEF ground state doublet by the Sm 
magnetic interaction/order. 
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FIG. 15. Parameters with best fit of equations (8) to the 
data. The fits are shown in Fig. 13 and 14. We have added 
the energy of the lowest CEF excitation determined by Cp{T): 
Al =20.0(2) meV for 2;=0.10 by Baker et al. [62] (diamond), 
and Al =22.0(2) meV for x—O.Od in this work (triangle). 



as a reference and proceeded as Baker et al. 
[62]. They reported three crystal electric field dou- 
blets for SmFeAsOo.9oFo.io with energies Ai^2=20.0(2), 
45(1) meV of the first and second excited levels, respec- 
tively. For .T=0.06 we find Ai,2=22.7(2), 52.7(9) meV. 
We plot Al as a function of doping together with A (from 
^SR) in Fig. 15. We find that A is systematically salO% 
larger than Ai. Part of the deviation should be due to 
the uncertainties of the analysis of both the specific heat 
and the fiSK data. In particular our model (8) neglects 
the second CEF excitation, one-phonon scattering, and 
the contribution of Sm-Sm magnetic interactions to the 
line width of the quasi-elastic CEF excitations. 



Superconductivity 

To investigate the magnetic penetration depth of 
SmFeAsOo.9oFo.io in the superconducting phase we ap- 
plied a transverse magnetic field ^qH=70 mT that causes 
muon spin precession. In general we expect the relax- 
ation above Tc to be solely due to the spin fluctuations 
described in Sec. . For a powder sample of an anisotropic 
type-II superconductor, the flux line lattice causes an ad- 
ditional Gaussian relaxation rate asc{T)=y^a{T)'^ — ct^^ 
where anm is determined above Tc [74]. age is a mea- 
sure of the superfluid density, or the inverse penetration 
depth (Tscocl/A^ocns/TO* [74]. Both effects contribute to 
the damping of the precession signal, which is modeled 

by 



A{t) = (ae" 



)e ^ cos{jf,Bioct- 



In conclusion, 



interpret A 



L.iA/B) 



the 



spin- 



_ (9) 

In the limit of fast spin fluctuations, the longitudinal and 
transverse relaxation are identical to Ai^(yi/s) [49]. Both 
can be measured simultaneously at the GPS instrument 
of the Paul-Scherrer Institute, which improves the ac- 
curacy of cr. For briefness we do not show \L,{A/B)iT) 
(they are identical to the rates measured in zero field, see 
Fig. 13). 

Bulk superconductivity is indicated by the reduction 
of the internal field Bioc and the increase of cFsc below 
Tc shown in Fig. 16. The saturation of (Jgc indicates 
a low density of states at the Fermi level characteristic 
for a nodeless gap. A gap with nodes, as e.g. for ex- 
tended s-wave symmetry would lead to a linear decrease 
of A at low temperatures [75]. asc{T) has been ana- 
lyzed using the temperature dependence of 1/ A^ assum- 
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ing a nodeless s-wave gap A(T) [76]. The best fit, shown 
in Fig. 16, yields A(T = 0)=10.0(5) meV, Tc=50.6(4), 
and crsc(0)=1.88(5) fis^^. For an anisotropic supercon- 
ductor the in-plane penetration depth can be estimated 
to Aa6(0)[nm]=248.85(crsc(0)[Ai s-i])-°-5=180(2) nm [74, 
77]. Vortex lattice disorder can artificially reduce Aab, 
therefore 180(2) nm is a lower limit for 180(2) nm. Nev- 
ertheless, this value is in agreement with the (more re- 
liable) results of Weyeneth et al. [78]. They studied 
polycrystalline SmFeAsOOssFo.is (nominal composition, 
Tc = 54 K from electrical resistivity) by TF-/xSR and 
found Xab ~ 200 nm by carefully measuring the field 
dependence of ■ This procedure allowed them to ex- 
clude vortex lattice disorder due to pinning. They were 
able to determine the sizes of both superconducting gaps 
A(0) = 13.8,5.3 meV. The resulting asc{T), shown as a 
dashed line in Fig. 16, is compatible with our data. How- 
ever, within the error of our data, a fit with two gaps is 
not possible. 
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FIG. 16. Weak transverse field /iSR of SmFeAsOo.9oFo.io: 
The Gaussian relaxation rate asc due to the flux line lattice, 
and the reduction of the local magnetic field Bioc due to the 
diamagnetic response below Tc indicate bulk superconductiv- 
ity. 



Summary 

In Tab. I we summarize the characteristic temperatures 
relevant for the structural and electronic phase transi- 
tions of SmFcAsOi-ajFa; determined by muon spin relax- 
ation, synchrotron X-ray diffraction, electrical resistivity, 
and specific heat measurements. 

General results The phase diagram is shown in 
Fig. 17. As a result of the detailed examination of the 
(broadened) phase transitions with /xSR and synchrotron 
X-ray diffraction the hierarchy of phase transition is 
not as clear as thermodynamic probes suggest (see e.g. 
Refs. 79 and 80). We observe long range magnetic order 
(Irmo) in clusters below Tsro~TP^^^>Ts already above 
the structural transition but the order parameter follows 



the expected temperature dependence only within the or- 
thorhombic phase {Tn<Ts). Both phase transitions are 
suppressed by doping and only the magnetic order, albeit 
now truly short ranged, survives as superconductivity ap- 
pears for a;>0.06. At higher doping levels we observe only 
superconductivity. The magnetic penetration depth of 
the superconducting material x=0.10 is consistent with 
a single nodeless s-wave gap. However, within the uncer- 
tainty of the analysis it is also possible that two gaps are 
present as suggested in several works [14, 78, 81, 82]. The 
magnetic transition of the Sm moments is almost unaf- 
fected by doping and is only reduced by wl K over the 
whole phase diagram. The spin fluctuations of the Sm 
are also present throughout the phase diagram and are 
almost unaffected by the disappearance of Fe magnetic 
order. 
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FIG. 17. The electronic phase diagram of SmFeAsOi-^Fa;. 
The extend and shape of the region with a mixed phase 
with superconductivity, Fe magnetic order, and non-zero or- 
thorhombicity is still unclear. In this work we only find 
a mixed phase of superconductivity and Fe magnetic short 
range order with almost identical transition temperatures that 
has a tetragonal crystal structure. This indicates that a;=0.06 
is close to a tetra critical point. Ref. 9 indicates that a re- 
gion with Fe magnetic long range order and superconductivity 
exists. Refs. 11, 68 indicate that superconductivity and Fe 
magnetic (long range) order occur also in the orthorhombic 
phase. All lines are guides to the eye. Dotted lines are extrap- 
olations of the phase boundaries indicated by solid lines. For 
a;=0.08, and 0.10, it is unclear whether Sm order occurs (see 
Sec. , p. 12), this is indicated by the open symbols and the 
dashed line. For 2;=0.06 the triangle that marks Tsro{~Tc) 
was drawn open for clarity. 

Disorder The structural phase transition is also not 
sharp, but occurs over a broad temperature interval, 
which may be caused by random or homogeneous stress, 
or order parameter fluctuations in a broad temperature 
interval. In this regard it is worth mentioning, that also 
the possible spin-nematic phase transition can be broad- 
ened by stress[83, 84]. The broad magnetic transition 
may also be due to homogeneous [44] or inhomogeneous 
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stress[83]. Vavilov and Chubukov [85] showed in a theo- 
retical study that disorder may have the same effect on 
the phase transitions as charge doping. In light of this 
result we cannot dismiss quenched disorder (impurities, 
lattice defects) as a second possible source of the broad 
phase transitions. 

Structural transition The tetragonal to orthorhombic 
phase transition does not occur for a;=0.06 where short 
range magnetic order occurs together with superconduc- 
tivity. This is supported by the absence of any extra 
anomaly of the temperature dependence of the specific 
heat Cp{T) (Fig. 11, p. 12). Also our X-ray diffraction 
study shows no indication for a structural phase tran- 
sition. It follows that the structural phase transition 
be absent for x>0.06. Due to the large FWHM of the 
(2,2,0)t Bragg peak, we are not able to resolve small or- 
thorhombic distortions (5<1T0~'^ (the discussion of the 
resolution limit on p. 21 also applies here). We note that 
small changes of the FWHM alone cannot be taken as 
proof for a structural phase transition, as this can lead 
to false phase transition temperatures (see p. 6). 

Mixed phase For a;=0.06 we find magnetic short range 
order below Tsro=33(l) K and superconductivity with 
rc=32.0(l) K. Whether the two orders coexist micro- 
scopically or phase separated cannot be decided at this 
point — the relevant evidence are: 

• bulk superconductivity (from x(T)[20]), and 

• bulk magnetic order (from ^SR, see pp. 10/). 

For microscopic coexistence both order parameters would 
be non-zero in the bulk and, possibly, we could observe 
a coupUng of order parameters if Tc<Tsro[^0, 56, 86, 87]. 
In case of phase separation the volume fractions should 
be coupled, i.e., one grows at the expense of the other[52- 
55, 88]. We also have to consider nanoscopic phase sep- 
aration that has been mentioned in recent ^SR works [7- 
9]. In this case even /iSR is not able to directly detect 
the phase separation. The available experimental evi- 
dence are compatible with both situations! Note, that 
concerning the question of nanoscopic phase separation 
vs. microscopic coexistence, our experimental results are 
consistent with the results in Ref. 9. We merely point 
out that this issue cannot be resolved using the available 
experimental data. 

Relevance for a possible QCP This ambiguous evi- 
dence leaves the question open whether a magnetic quan- 
tum critical point could occur in the phase diagram of the 
Sm system, as it is not clear whether a second or first or- 
der phase transition occurs between magnetic order and 
superconductivity as a function of doping. 

Electrical resistivity & phase transitions The rela- 
tionship between the electrical resistivity and the struc- 
tural, and magnetic phase transitions is illustrated in 
Fig. 18. Apparently, no simple connection exists be- 
tween the magnetic, and structural phase transitions and 



temperature dependence of the electrical resistivity. The 
maximum of the resistivity at T=T^^^ does not co- 
incide with any phase transition studied in this work 
since it occurs at somewhat higher temperatures. The 
maximum marks the onset of a reduction of the critical 
fluctuations prior to the structural/magnetic transitions, 
and thus indicates the development of short range struc- 
tural/magnetic order. The inflection point agrees with 
the measured Tm for a;<0.04 within errors. The upon 
coohng reduced decrease of p{T) signals an effective de- 
pletion of electronic states at the Fermi level which is ex- 
pected for a true long-range ordered spin density wave. 
The mismatch between Tf and T^r could be related to 
the width of the transition. 

CeFeAsOi_;rF:r 

In this section we present the results of our electrical 
resistivity (Sec. ), X-ray diffraction (Sec. ), muon spin re- 
laxation measurements (Sec. ), Mossbauer spectroscopy 
(Sec. ), and magnetic susceptibility measurements of the 
materials that is both superconducting and magnetically 
ordered (Sec. ). Our experimental data show that this 
system is similar to the Sm system and in most parts 
we applied the same analysis strategies and draw simi- 
lar conclusions as in the preceding section. In Sec. we 
summarize our main results and present the electronic 
phase diagram of CeFeAsOi-^Fa:. All phase transition 
temperatures are summarized in Tab. III. 

Electrical resistivity measurements 

CeFeAsO shows, similar to SmFeAsO, a characteris- 
tic cusp in the temperature dependence of the electrical 
resistivity p{T) at T^^^^^lbO K followed by an inflec- 
tion point at 7;''=135(2) K in Fig. 19. For the doped 
compounds, this anomaly is not observed, presumably 
due to the higher doping level compared to the non- 
superconducting Sm-based compounds. The suppression 
of the cusp is correlated with the reduction of Tg and 
the low temperature limit of the orthorhombic distortion 
S{T—>-0) to values smaller than those we presented for 
SmFeAsOo.96Fo.o4- Instead of the cusp, we observe a re- 
sistivity minimum at T«70 K. Later we will show that 
Ts is above, and T/v close to the minimum of p(T) (see 
Sec. ). For a:=0.048(3), p{T) is similar to a;=0.042(2). 
However, the minimum is followed by a cusp at Tsa26 K. 
At lower temperatures the resistivity drops only almost 
to zero which indicates spurious superconducting grains 
and not bulk superconductivity. 

We observe bulk superconductivity first for a;=0. 063(2) 
below rc=29.5 K. At higher temperatures p{T) shows a 
cusp at rw35 K and a minimum at Tk70 K much like 
for a;=0.048(3). Above rw240 K p{T) becomes linear in 
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FIG. 18. A comparison of the X-ray diffraction patterns, orthorhombic distortion 5 and the FWHM of the (2,2,0)t Bragg peak, 
the electrical resistivity p and its derivative, and the muon spin precession frequency transverse relaxation rate Xt/{2tv) 
and magnetic volume fraction Vmag of the Sm samples with a structural phase transition. The vertical dashed lines indicate 
the phase transition temperatures. Where applicable arrows indicate the maximum (T^^^) and the inflection point {T[) of 
p{T). This comparison shows that the phase transitions cause the anomalies of the resistivity, but no simple relationship exists 
between the phase transition temperatures and the shape, and position of the anomalies. Except the color coding of the density 
plots, the j/-axes scaling is identical for the different doping levels. 



T. For a;=0.145(20), and 0.150(20), we find supercon- 
ductivity below rc=43.8, 43.4 K, respectively. At higher 
temperatures the resistivity shows neither a cusp nor a 
minimum but follows p(r)cxT^ for both materials. 

In the superconducting regime, the resistivity of 
CeFeAsOi-xF^: shows similarities to LaFeAsOi-xFa^, 
and SmFeAs0i-2:F2:. In the optimally doped and over- 
doped regime {x ^ 0.1) LaFeAs0i_a;F2, also shows 
linear p{T) at high temperatures and behavior at 



low temperatures [22]. In the underdoped regime of 
LaFeAs0i_2;Fa; (a;<0.1), p{T) is also proportional to 
temperature above Tw200 K [22]. At lower tempera- 
tures, p{T) drops below the (linear) extrapolation of the 
high-tempearture linear behavior. This resistivity drop, 
accompanied by a large Nernst effect anomaly, was inter- 
preted as the signature of slow spin density wave fluctu- 
ations [22, 89]. The x=0. 063(2) doped Ce sample shows 
the same feature in the resistivity which suggests a sim- 
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TABLE I. Phase transition temperatures for tlie structural transition {Ts, see Sec. , p. 3), superconducting phase transition (Tc 
obtained by /iSR and electrical resistivity p{T<Tc)~0, see Sec. , p. 15 and , p. 3), onset of magnetic order (Tgro), bulk magnetic 
long range order (TV), saturation of the magnetic volume fraction at r=rioo (Tjv, Tsro, Tioo obtained by pSR, see Sec. , p. 7). 
T^™ is the magnetic ordering transition temperature of the Sm moments obtained by specific heat measurements (cp), from 
the /iSR transverse relaxation rate (At), or the fiSK transverse relaxation rate (Al), see Sec. , p. 12. All temperatures are 
given in K. *The Sm transition temperatures determined from \l{T) should not be considered as proof for a phase transition 
(see Sec. , p. 12 for details). 




temperature T (K) temperature T (K) 

FIG. 19. Normalized resistivity as a function of tempera- 
ture of CeFeAsOi-iFi. (Blue and green) Red curves refer 
to (non-)superconducting compounds, right: b) the blue line 
is a linear fit above 200 K, below r=200 K resistivity dif- 
fers from linearity similar to underdoped LaFeAsOi-a:Fi:. c) 
Quadratic temperature dependence of p{T) below r~200 K in 
resemblance of overdoped LaFeAsOi-a:Fi. Curves are shifted 
vertically for clarity. 



ilar role of magnetic fluctuations. In Sec. we will show 
that, contrary to the La system, these fluctuations slow 
down and eventually lead to magnetic order. 

In the discussion in Sec. we will compare the struc- 
tural and magnetic phase transition as found by X-ray 
diffraction, /xSR and Mossbauer spectroscopy to the fea- 
tures in the temperature dependence of the resistivity. 

Susceptibility measurements 

The temperature dependence of the magnetic suscep- 
tibility for a;=0.063(2) in magnetic flclds of 0.805(3) mT, 
and 0.112(3) mT cooled in flcld (fc) and in zero flcld 



(zfc) is shown in Fig. 20. Both fc and zfc measurements 
show strong diamagnetic signals in the superconducting 
state. However, neither reaches Attx{T)=—\ down to 
T«2 K. The sizable fc and zfc diamagnetic susceptibil- 
ities indicate bulk superconductivity. Later we present 
X-ray diffraction (via the jump thermal expansion coeffi- 
cient) data that also support bulk superconductivity, and 
/iSR data that is consistent with nearly 100% supercon- 
ductivity for Tc>r>20 K. The granular, i.e., polycrys- 
tallinc structure of the samples may cause the reduced 
shielding signal: A'kx{T)=—1 may only be observed if 
the weak Josephson-junctions between the grains provide 
superconducting paths for the Meissner shielding current 
across the whole sample, e.g., see Ref. 90. Hence, we con- 
clude that the smallest applied field of about 0.1 mT, pos- 
sible with our experimental setup already breaks some of 
the weak links. However, the data clearly allow extract- 
ing Tc- Note, that the onset of Ce magnetic order at 
T^'^=2.7(l) K is visible in the data by the slight upturn 
and cusp, too. This is consistent with the measurements 
in Rcfs. 7 and 8, and our /xSR measurements (see below). 
Also 47rx(T)>— 1 is in agreement with results of Sanna 
et al. [7] and Shiroka et al. [8] Based on ^SR results, they 
also concluded bulk superconductivity in similarly doped 
materials with a comparable macroscopic susceptibility. 
Important is also the almost constant susceptibility for 
T<20 K, it shows that the superconductivity is not, and 
not even in parts, quenched by the low temperature Fe 
magnetic order present in this sample. 



X-ray diffraction measurements 

As for the Sm system, we study the structural phase 
transition of the Ce system using synchrotron X-ray 
diffraction. We find that it is qualitatively identical in 
both systems. The phase transition is also broad and 
preceded by lattice fluctuations. Here we only summa- 
rize the features. Please see Sec. , p. 3 for a detailed 
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discussion of our analysis strategy. 




temperature T (K) 



FIG. 20. The temperature dependence of the magnetic sus- 
ceptibility 47rx(T) of CeFeAsOi-:,F:j, with a;=0.063(2). Open 
symbols denote measurements upon heating after the sample 
was cooled in zero magnetic field (zfc) and closed symbols 
after cooling in the applied magnetic field (fc). 

Representative diffraction patterns of CeFeAsO, and 
the orthorhombic distortion are shown in Fig. 21. We 
determined the lattice constants a, 6, and c in Fig. 22 
from the positions of the tetragonal (2,2,Q)t and (0,0, 6)t, 
and the orthorhombic (4,0,0)o and (0,4,0)o Bragg peaks. 
The temperature dependence the orthorhombic distor- 
tion in Fig. 21, the lattice constants in Fig. 22. and 
the FWHM of the tetragonal peak above the transi- 
tion in Fig. 23 follow the same trends we reported for 
the Sm system. The tetragonal peak also broadens sig- 
nificantly when we approach the structural transition 
from high temperatures. We find rs=148.9(-h2-8) K, 
97.0(-hl-7) K, 74.7(+l-5) K for x=0, 0.042(2), and 
0.048(3), respectively (the positive error is due to the 
temperature step of the measurement and the nega- 
tive error is the temperature range for which 5{T) 
is convex). The splitting is apparent in the raw 
data shown in Fig. 21. It persists only for the non- 
superconducting compounds with a;=0, 0.042(2), and 
0.048(3), and is absent for the superconducting com- 
pounds with a;=0.063(2), and 0.150(20). 

At first glance, the structural transition tempera- 
ture Ts and the low temperature orthorhombic distor- 
tion (5(T— )-0) are strongly reduced in the two doped 
CeFeAsOi-xF^^ specimens as shown in Fig. 21. However, 
Fig. 24 indicates that the reduction of ^(T— >-0) is propor- 
tional to Ts, i.e., 5{T^Q,Ts)(xTs and that 5{T^Q,Ts) 
is quantitatively identical for the Sm and the Ce system. 
In Fig. 24 we also show the data obtained by Margadonna 
et al. [11] for the Sm system, and Zhao et al. [4] for the Ce 
system. For the non-superconducting Sm and Ce materi- 
als 5(T— >-0, Ts) follows the same linear dependence on Ts- 
For the two superconducting specimens studied by Mar- 



gadonna et al. [11] the orthorhombic distortion is higher 
than expected from the linear 5{T^Q,Ts) at lower dop- 
ing levels, whereas the superconducting specimen stud- 
ied by Zhao et al. [4] follows the linear behavior. At a 
quantum critical point the order parameter and the tran- 
sition temperature should simultaneously become zero. 
However, a linear extrapolation of (5(T— >-0, Ts), shown in 
Fig. 24 indicates that transition temperature and order 
parameter do not simultaneously become zero (exclud- 
ing the two outlying data points taken from Margadonna 
et al. [11]). On the other hand, the two outlying data 
points suggest that the decrease of the orthorhombic dis- 
tortion slows down, while Ts rapidly vanishes in the su- 
perconducting regime. Also in this case, both parame- 
ters do not simultaneously vanish. However, a structural 
QCP (or the related nematic QCP [91]) may be recov- 
ered by a "back bending" of the structural phase transi- 
tion line, which has been reported for doped BaFe2As2 
[10, 38, 92]. 

The analysis of the orthorhombic distortion with the 
general power law i5(T)cx(l — (T/Ts)")^ is, just as for 
the Sm system, not feasible due to the broadened transi- 
tion. But we can compare the temperature dependence 
of the normalized orthorhombic distortion with 2D-Ising 
and 3D-Ising order parameters (see Fig. 25). Qualita- 
tively, this indicates that the order parameter univer- 
sality class changes from 2D-Ising (x=0) to 3D-Ising 
(x=0. 048(3)) via an intermediate temperature depen- 
dence (x=0. 042(2)) similar to the Sm system (see Fig. 5). 

For a;=0. 063(2), and 0.150(20) we observe an upturn of 
the a lattice constant below Tk,Tc in Fig. 26 (c/. Sec. ). 
We can exclude a magnetic origin of the upturn (see 
next sections): for x=0. 063(2) a sizable magnetic vol- 
ume fraction develops only at lower temperatures, and 
x=0. 150(20) is non-magnetic. The expansion of the a 
axis shown in Fig. 26 is most likely due to supercon- 
ductivity, and suggests AaaOcdTc/ dpa < (the effect of 
superconductivity on the c lattice constant could not be 
resolved). More importantly, this finding for x^O. 063(2) 
supports that a substantial part of the sample becomes 
superconducting. 

The FWHM (fuU-width-at-half-maximum) of the 
(2,2,0)t peak in Fig. 23 for a;=0. 063(2) strongly increases 
at low temperatures. It reaches a broad maximum at 
TwlOO K. This coincides with an inflection point in the 
electrical resistivity in Fig. 19. The minimum of p{T) 
at rsa70 K roughly coincides with the minimum of the 
FWHM of the (2,2,0)t Bragg peak. For a;=:0. 150(20) the 
FWHM shows no sharp anomalies, and also p{T) shows 
no pronounced features. Similar to the Sm system, the 
absolute value of the FWHM is reduced by «50% com- 
pared to a::=0. 063(2). Contrary to the Sm system, we 
observe no anomalies at Tc for both a;=0. 063(2), and 
0.150(20) which speaks against suppression of lattice fluc- 
tuations by superconductivity. 

The question arises whether the broadening of the 
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FIG. 21. Left: representative X-ray diffracted patterns of the (2,2,0)t tetragonal and the (4,0,0)o and (0,4,0)o orthorhombic 
Bragg peaks at different temperatures for CeFeAsO. Right: The orthorhombic distortion as a function of temperature and 
doping level of CeFeAsO i- a; F^:. Gray data points indicate the temperature range for which 5{T) is convex (see Sec. , p. 3 for 
a detailed explanation). 
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FIG. 22. The lattice constants as a function of temperature 
of CeFeAsOi-iFi,. The black circles refer to a and b lattice 
parameters and the red triangles to c. Below the tetragonal 
to orthorhombic transition, a^b are divided by for com- 
parison. 
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FIG. 23. FWHM of the (2,2,0)t (left) above the structural 
transition, and (right) in the superconducting superconduct- 
ing regime of CeFeAsOi-iFa, arrows indicate the structural 
and superconducting transition temperatures Ts and Tc and 
the minimum of the resistivity (T^^^) for a:=0. 063(2). 



tetragonal Bragg peak for a;=0. 063(2), and 0.150(20) in 
Fig. 23 could be an unresolved orthorhombic splitting. 
The minimum and the following increase of a;=0. 063(2) 
for r<70 K could indicate an unresolved splitting (see 
Fig. 23). To answer this question, wc discuss the resolu- 
tion limit of our X-ray diffraction experiment below. It 
is mainly limited by the intrinsic width of the (2,2,0)t 
Bragg peak. For T=70 K, the (2,2,0)t Bragg peak is at 
26'=45.02° with a FWHM of «0.1088°. The angle resolu- 
tion of our experiment is «0.013° and the error of each in- 
tensity value is «3%. To estimate the resolution limit we 
simulate a series of diffraction patterns of a hypothethi- 
cal orthorhombic phase with two Gaussian Bragg peaks 
((0,4,0)o and (4,0,0)o) at 26'±=45.02°±0.5A with the 



FWHM of the tetragonal peak (0.1088°) and the above 
angle resolution. We fit each simulated pattern with a 
single Gaussian Bragg peak. It turns out, that for an 
orthorhombic distortion 5=0.40-10-3 (A=0.0381°) the 
simulated diffraction pattern is indistinguishable from a 
single Bragg peak with a FWHM of «0.12°, which corre- 
sponds to the experimental FWHM at TwlO K. There- 
fore, S<QA0-1Q^^ is the upper limit of the orthorhombic 
distortion that is still consistent with the experimental 
data. To unambiguously resolve the orthorhombic dis- 
tortion, its error should be smaller than the orthorhom- 
bic distortion itself. With the given FWHM and exper- 
imental resolution, this is satisfied for an orthorhombic 
distortion larger than (Jwl-lO^'^. 
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FIG. 24. The low temperature limit of the orthorhombic dis- 
tortion 5(r— >-0) as a function of the structural transition tem- 
perature Ts for the Sm and the Ce system. The dotted line 
is a linear fit and the fit parameters are given in the figure. 
For comparison, we also show the data taken from the work 
of Zhao et al. [4], and Margadonna et al. [11]. The mark- 
ing on the a;-axis indicates TJ! for which 5{T ^ 0,T^) = 0, 
and the bars indicate the error of Tj (one, and two standard 
deviations) . 
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FIG. 25. The orthorhombic distortion 6{T)/S{T=0.6Ts) nor- 
malized to the value at T=0.6Ts as a function of reduced 
temperature T/Ts for CeFeAsOi_a;Fi,. The Sm and the Ce 
compounds show similar doping and temperature dependen- 
cies (compare with Fig. 2, p. 4). To illustrate the change 
of the universality class of the order parameter we also show 
the expected temperature dependence for an order parame- 
ter of the 2D-Ising (oc(l - T/Tsf''^^^), and the 3D-Ising type 
(oc(l — T/TsY'''^^^). Due to the large asymmetric error of 
Ts we could not determine the universality classes unambigu- 
ously. Within the uncertainty of the analysis a;=0.00 belongs 
to the 2D-Ising universality class, whereas a;=0. 048(3) belongs 
to the 3D-Ising class and 2;=0. 042(2) lies between the two. 
The change of dimensionality is more pronounced in the Ce 
than in the Sm system. This is possibly due to the lower 
transition temperatures of the Ce materials. 



FIG. 26. The low temperature upturn of the a lattice constant 
below T«Tc of CeFeAsOi-xF;i; with a;=0. 063(2) (top), and 
a;=0. 150(20) (bottom). This is similar to the upturn of the 
a and c lattice constants of SmFeAsOi-^F^ with a::=0. 048(3) 
(see Fig. 6, p. 5). It is due to the pressure dependence of Tc. 
The change of slope from negative {T<Tc) to positive {T>Tc) 
suggests that Tc would decrease under in-plane pressure. We 
were not able to resolve the effect of superconductivity on 
the c lattice constant with our X-ray diffraction experiments. 
Because the dependence of Tc on hydrostatic pressure depends 
both on in-plane (along a) and out-of-plane (along c) pressure 
we can draw no conclusion as to whether Tc should increase 
or decrease under hydrostatic pressure. 



Several studies of ferropnictides suggest [79, 93, 94] that 
the orthorhombic distortion is proportional to the or- 
dered Fc magnetic moment. A linear extrapolation using 
the orthorhombic distortion and the Mossbauer Fe hyper- 
fine field, which is directly proportional to the ordered Fe 
moment (see Sec. ) suggests an orthorhombic distortion 
of «0. 5(1)10^'^ for .T=0. 063(2) at lowest temperatures. 
This value is still compatible with the data but our ex- 
periment would not be able to resolve this distortion. 



In summary, the resolution limit of our diffraction 
experiment is S^l-lO^^ and the largest orthorhom- 
bic distortion that is still compatible with the data 
is (5sa0.4-10~'^. Considering the above estimates of 
(5(r — > 0) we cannot exclude a structural transition for 
x=0. 063(2). However, in that case, 6{T 0) should be 
smaller than w0.4-10~"^. 
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FIG. 27. Zero and longitudinal field muon spin relaxation time spectra of CeFeAsOi-ajFj, and (bottom right) a weak longitudinal 
field decoupling time spectrum for 2;=0. 063(2). Note the different time scales for the various plots. 
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FIG. 28. The muon spin precession frequency is propor- 
tional to the AFM order parameter and indicates long range 
magnetic order. 



Muon spin relaxation measurements 

Fe magnetic order 

Zero field fiSR data for x<0. 063(2) are best described 
by a model for commensurate magnetic order; 

Ait) =Vmag [(1 - ai)a2(^e-^^' cosLot + 

+ (l-a2)e-^^i (10) 
+ (l-V;„ag)G(t,a)e-^-. 



The first term models the commensurate Fe magnetic or- 
der. The second term is non-zero (a2<l) for the doped 
materials only below sa4 K where Ce magnetic order 
causes an additional fast relaxation with rate \ce due 
to magnetic fluctuations evident from the reduction of 
P{t) below 1/3. The last term models the paramag- 
netic phase with weak Gaussian-Kubo-Toyabe relaxation 
G{t,a), and magnetic fluctuations that cause relaxation 
with rate A„m [49]. The fluctuations close to the mag- 
netic transition cause part of the signal to relax during 
the dead time of the detectors, this is modeled by the 
factor ai that is only non-zero and of the order of a few 
percent close to the magnetic transition. 

When we cool the non-superconducting compounds 
with a;<0. 048(3) the magnetic order manifests in the ^SR 
data in two steps: First magnetic short range order devel- 
ops at r^™=145(5) K, 75(3) K, 49(1) K and the magnet- 
ically ordered phase volume Vmag starts to increase grad- 
ually until it reaches 100% at rioo=133(2) K, 40(2) K, 
30.0(5) K for a;=0,0.042(2), and 0.048(3), respectively. 
This can been seen by the absence of an oscillation and 
the fast exponential relaxation with rate Xt ^ 20 MHz 
of the muon spin polarization A{t) for t < 0.5 /iS that 
reaches a signal fraction of « 2/3 when the material is 
fully ordered (see Fig. 27). This indicates that magneti- 
cally ordered clusters with short coherence lengths form 
throughout the material below Tsro that saturate at Tioo. 
We observe coherent long range magnetic order with 
a well defined order parameter, the muon spin preccs- 
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FIG. 29. The magnetic volume fraction Vmag with fits according to power law distributed transition temperatures. A small 
but non-zero Vmag indicates magnetic short range order. 



sion frequency only close to Tioo with TAr=137(2) K, 
37(1) K, 32.0(5) K for x=0, 0.042(2), and 0.048(3), re- 
spectively. This resembles the broadened magnetic phase 
transition of SmFeAsOi-^^F^ that we reported in Sec. . 
However, in the Ce system, we find no coherent preces- 
sion signal above Tioo which suggests that the magneti- 
cally ordered clusters arc too small or too disordered for 
long range magnetic order. f^{T) shown in Fig. 28 is 
approximately four times smaller in the doped materi- 
als, similar to the observations [5, 9] in SmFcAsOi-^rFx 
in Sec. . This is different than in LaFcAsOi-jjFa: where 
fix is only reduced by 10% before superconductivity is 
induced by doping [6]. 

The width of the transition AT—Tsro-Tioo in the un- 
doped material is 6(2) K and increases drastically in the 
doped materials up to 32(2) K for a;=0.042(2). This 
broadening is not symmetric, VmagiT) follows the same 
shape that we reported for the Sm system (see Fig. 29, 
and Sec. ). To quantify the broadening we model the 
VmagiT) with the distribution of ordering temperatures 
(5) and (3). Best fits are shown in Fig. 29. Opti- 
mized parameters are ro=90(8) K, 53(4) K, 37(6) K, and 
cr=36(8) K, 16(4) K, 25(6) K for a;=0.042(2), 0.048(3), 
and 0.063(2), respectively (a indicates the width of the 
transition, VmagiTo — (t)=0.5). The exponent 6=1.8(7) is 
independent of doping. For the Sm system we obtain a 
larger exponent 6=4.0(5), interestingly this value is close 
to what Bcrchc et al. [61] determined for the McCoy- Wu 
model: Vynag{T)(xT'^-'^ . Within this work, we did not de- 
termine Vmag{T) accuratc enough close to the onset of 
magnetic order to accurately determine the exponent 6. 
In particular, we cannot exclude the possibility of expo- 
nential decrease of Vmag{T) towards high temperatures, 
which would point to statistically rare unperturbed do- 
mains (for a review sec Ref. 95). On the other hand, the 
exponent determined by our analysis of Vmag{T) with a 
power law may be related to the type of disorder or the 
underlying magnetic interactions. 

We also find a broad magnetic transition for 



a;=0. 063(2), as can be seen by a non-zero Vmag below 
Tsro = 28(1) K in Fig. 29. But, corresponding com- 
pound of the Sm system, we detect no coherent muon 
spin precession which indicates that this is short range 
order. To exclude that the fast relaxation for t < 0.5 ^s 
in zero magnetic field is due to magnetic fiuctuations we 
measured in an applied longitudinal magnetic field (LF). 
It is well known, that a small LF will not have a signifi- 
cant cff'cct on the muon spin relaxation rate if it is solely 
due to magnetic fluctuations. Hence, the decoupling in 
LF seen in Fig. 27 is a clear indication for predominantly 
static magnetic order that develops in the material below 
Tsro = 28(1) K. 

For x=0. 145(20), we find, for a paramagnetic material, 
an unexpectedly fast relaxation of the muon spin polar- 
ization shown in Fig. 27 — but it is wlOO times smaller 
than for any of the magnetically ordered materials. This 
relaxation is temperature independent between 50 K and 
1.6 K and can be decoupled in a longitudinal magnetic 
field of 70 mT, see Fig. 27. This indicates that it is caused 
by dilute magnetic impurities. This is different from the 
relaxation in the superconducting Sm compounds where 
magnetic fluctuations of the Sm 4/ moments dominate 
(see Fig. 13, p. 14). This is consistent with the w5 wt.% 
of Fe2As impurities detected by Mossbaucr spectroscopy 
(see next section) 

For x=0. 150(20) we find a weak relaxation of Gauss- 
Kubo-Toyabe[49] form, typical for a non- magnetic ma- 
terial. Only at lowest temperatures we observe an ad- 
ditional exponential relaxation, most likely due to the 
slowing down of Ce magnetic moment fluctuations, see 
Fig. 27. 

Ce magnetic order 

In the undoped material the onset of long range Ce or- 
der is marked by a cusp in /^(T). This is also visible for 
a;=0.042(2), and 0.048(3) but less pronounced. This sug- 
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gests long range Ce magnetic order below T^'^=4.4(3), 
3.5(10), 3.5(10) K for a;=0, 0.042(2), and 0.048(3), re- 
spectively. The larger error for the doped material stems 
from the large relaxation rate Xt which increases the er- 
ror of and masks the onset of the Ce order. It is 
accompanied by an extra signal fraction that is fully dy- 
namic, i.e., due to rapid magnetic fluctuations (a2<l in 
(10)). It can be seen in the time spectra in Fig. 27 by 
the lowering of the polarization for t>0.25 us below 1/3 
expected for static magnetic order. 

For higher F content, long range Fe order is suppressed 
and the Ce order is weakened. For x=0. 063(2), the Ce or- 
der is marked by aforementioned dynamic signal fraction 
at r^'^=2.5(10) K which is accompanied by a cusp in the 
static relaxation rate At, i.e., the width of the local mag- 
netic field distribution. For x^O. 145(20) T^'==1.0(2) K 
is further suppressed. Above T=l K impurities dom- 
inate the muon spin relaxation (see previous section). 
Below T = 1 K the relaxation rate increases by a fac- 
tor of 3, marking the onset of Ce magnetic order. For 
x=0. 150(20) the relaxation rate is further suppressed to 
0.01 /is^^ which is of the order of the relaxation rate 
caused by randomly oriented nuclear magnetic dipoles. 
This rate increases to 0.03 fj,s~^ at T«2 K, which is too 
small to be due to magnetic order. This means that we 
do not observe static magnetic order of the Ce subsystem 
for the a;=0. 150(20) sample down to the lowest temper- 
ature of T=1.7 K measured in this experiment. 

In summary, we observe no long range Ce magnetic or- 
der in the absence of Fe magnetic long range order. The 
magnetic ordering of the Ce subsystem may stem from 
RKKY interaction between the localized Ce-4/ moments 
mediated by the Fc conduction band electrons. Akbari 
et al. [96] showed that the RKKY interaction is sizable in 
the fcrropnictidcs and that doping, Fermi surface topol- 
ogy, and the SDW ordering of the Fe moments influence 
the RKKY interaction of local moments. Additionally, 
disorder induced smearing of the conduction bands may 
also reduce the RKKY interaction strength and lead to 
a destruction of long range correlations and a reduction 
of the ordering temperature. 



Superconductivity 

As for the Sm system we use transverse magnetic field 
/iSR to determine the magnetic penetration depth. To 
account for the magnetic order for x^O. 063(2), we intro- 
duce a corresponding signal to the fit function of the fiSK 
time spectra: 
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FIG. 30. The temperature dependence of the parameters of 
the SC phase, top the volume fraction of the short range order 
of a;=0. 063(2), middle: the internal field that is reduced by 
the diamagnetic response, bottom: the SC order parameter 
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The first term models the magnetically ordered phase vol- 
ume, and the second term describes the (purely) super- 
conducting phase volume (c/. Sec. ) The close agreement 
between Vmag (T) obtained from the weak transverse and 
zero field /iSR in Fig. 30 indicate that the description of 
the /iSR spectra is consistent. asc{T)—\/cr{T)'^ — (^nm 
a measure of the superfluid density, or the inverse pen- 
etration depth cTscOcl/A^ocris/TO* [74]. The relaxation 
rate At/ accounts for relaxation due to static and dy- 
namic magnetism in addition to the relaxation due to 
the fiux-line lattice. 

For a;=0. 063(2), the diamagnetic reduction of the local 
magnetic field, and the additional Gaussian relaxation 
shown in the middle and bottom panel of Fig. 30 indicate 
bulk superconductivity with Tc=29.8(2). Close to Tc, the 
magnetic volume fraction Vmag{T) docs not grow larger 
than «10% for r>17 K, see top panel of Fig. 30. 

For a;=0. 145(20), static magnetism, due to dilute mag- 
netic impurities is already present above for T>Tc and 
in zero field the relaxation rate does not change between 
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50 K and 1.6 K (see Fig. 27, p. 23). In the applied 
transverse field ^qH=70 mT, Xtf{T) increases linearly 
from 0.7 ^ls~^ at r=50 K to 1 /^s^^ at r=3 K. For 
a;=0. 150(20), the transverse relaxation can be described 
fully by setting Knog=0. This indicates a homogeneous 
superconducting phase with dilute magnetic impurities. 
The bulk superconductivity is confirmed by the reduction 
of Bscc and the increase of CTsc for T < Tc=39.G(9). 

For a;=0. 150(20), in addition to the Gaussian relax- 
ation we observe a temperature independent slow relax- 
ation with At/(T)=0.05(1) iJ,s~^ due to fluctuations (see 
previous sections). They are enhanced by applied mag- 
netic fields >0.1 T and also present above Tc- The reduc- 
tion of Bsc and the increase of (Tsc at Tc=39.0(9) indicate 
bulk superconductivity, as above. 

The saturation of age for a;=0.2, and 0.150(20) indi- 
cates a low density of state at the Fermi level character- 
istic for a nodeless gap. A gap with nodes, as e.g. for ex- 
tended s-wave symmetry would lead to a linear decrease 
of A at low temperatures [75]. For all three supercon- 
ducting materials (Tscocl/A^ocris has been analyzed using 
the temperature dependence of 1/ A^ assuming a nodeless 
s-wave gap A(T) [76]. Fits are shown in Fig. 30. Best 
agreement has been found with the fit parameters shown 
in Tab. II along with the in-plane magnetic penetration 
depth Aafc(O). This closely agrees with the roughly linear 
scaling of Tc with the superfluid density ngOcl/A^ ob- 
served in many cuprate and ferropnictide superconduc- 
tors (see Fig. 35). Since we performed our measurements 
on powder samples in which vortex lattice disorder can 
result in an additional damping of the fiSK signal, which 
in principal cannot be disentangled from asc, the values 
listed in Tab. II have to be regarded as lower limits of 
the penetration depth. 



X 




A(0)(meV) 


2A(0) 
ksTa 




A„5(0)(nm) 


0.063(2) 
0.145(20) 
0.150(20) 


29.8(2) 
39.0(9) 
38.7(1) 


4.5(1) 
5.2(2) 
5.06(2) 


3.5(1) 
3.1(1) 
3.04(1) 


0.54(1) 
1.94(4) 
1.93(1) 


339(5) 
179(5) 
179(1) 



TABLE II. Best fit parameters of a^dT) oc 1/A(r)^ for a 
nodeless s-wave gap. 



Mossbauer spectroscopy measurements 

In the following, we present our ^^Fe Mossbauer spec- 
troscopy experiments. The results corroborate the re- 
sults of the /xSR experiments presented above. Typical 
Mossbauer absorption spectra are shown in Fig. 31(a). 
We recorded the spectra using a ^''Co/Rh gamma radi- 
ation source (emission line half-width-at-half-maximum 
r=0. 130(2) mm/s) and a silicon drift diode detector. We 
calibrated the spectrometer to a-Fe at room temperature. 



At room temperature, all spectra of CeFeAs0i_2;Fa; 
show a single sharp Lorentzian absorption line with 
a HWHM r=0.15(l) mm/s and an isomer shift of 
(5^^=0.45(1) mm/s with respect to a-Fe. This is con- 
sistent with the results of other Mossbauer studies of 
i?FeAsOi_^F^[23, 60]. 

Magnetically ordered Fe bearing impurity phases can 
be detected with Mossbauer spectroscopy if the con- 
centration of the impurity is larger than wl atom%. 
Only the spectra of superconducting material with 
a;=0. 145(20) show absorption lines in addition to the 
main line. The additional absorption lines are consistent 
with the Fe hyperfine field[97] of Fe2As and their spectral 
weights are equivalent to «5 wt.% of Fe2As impurities. 
The spectra of x=0, 0.048(3), and 0.063(2) show no addi- 
tional absorption lines at room temperature, which indi- 
cates that these samples do not contain any iron-bearing 
impurities. 

The spectra of a;=0.00 show a distinct broadening of 
the singlet pattern below 135 K that develops into a re- 
solved sextet at low temperatures due to the Fe magnetic 
order consistent with other studies of i?FeAsO[23, 60]. 
We analyzed the spectra of a;=0.00 like in Ref. 23 by 
diagonalizing the hyperfine Hamiltonian including elec- 
tric quadrupole and static magnetic hyperfine interac- 
tion. To account for the small broadening we used two 
sextets to describe the absorption spectra. The tempera- 
ture dependency of the magnetic hyperfine field Bhf (see 
Fig. 31(b)) and the low temperature values of the iso- 
mer shift (5^^=0.58(1) mm/s, the quadrupole shift e=- 
0.04(1) mm/s and the Fe hyperfine field Bhf=5.6(l)T at 
T=12 K are consistent with the results of McGuire et al. 
[60]. In addition to the hyperfine field of the Fe order 
they found that the magnetic order of the Ce moments 
increases the measured Fe hyperfine field slightly below 
«5 K due to a small transferred hyperfine field. 

At low temperatures, the absorption lines of the doped 
compounds with a;=0. 048(3), and 0.063(2) broaden sig- 
nificantly at low temperatures. Down to w4 K we observe 
no resolved sextet as for CeFeAsO. The broad absorption 
lines indicate Fe magnetic order with a small ordered mo- 
ment and a distribution of ordered moments. 

We analyzed the spectra of a;=0. 048(3) with a Gaus- 
sian distribution of hyperfine fields i?hf with width a and 
mean (i3hf), shown in Fig. 31. It turns out that the ratio 
(-Bhf )/o'= 1.41 is temperature independent. 

For a;=0. 063(2) this analysis is not possible because in 
this case the fit parameters a, and (Bhf) are too strongly 
correlated. Instead we fit a single Lorentzian line to the 
absorption spectra and extract the magnetic contribution 
,5HWHM(T)=0.5(r(T)-r(300 K)) in Fig. 31 to the total 
line width r(r) by subtracting the line width at room 
temperature r(300 K). For a quantitative comparison of 
(5HWHM with the hyperfine field of a;=0, and 0.048(3) 
we additionally analyzed the spectra of a;=0. 048(3) with 
a single Lorentzian line. It turns out, that for x=0. 048(3) 
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FIG. 31. (a): Typical Mossbauer spectra. The full lines are fits to the data and corresponding subspectra (see text). (6): 
Mossbauer Fe magnetic hyperfine field for a;=0, and 0.048(3) and the magnetic contribution 5HWHM to the absorption line 
width for x^O. 063(2). (c): /iSR transverse relaxation rate At (2;=0. 048(3), 0.063(2) are emphasized for clarity). Lines are 
guides to the eye only. Note the consistency of data produced with /iSR and Mossbauer spectroscopy in (6) and (c). 



we recover the ratio (Bhf)/(5HWHM=14.7 T/(mm/s), 
which is specific to the ^^Fc nucleus, from our fits for 
all temperature. This proves the validity of our analysis 
and allows us to compare (5HWHM of a;=0. 063(2) with 
the Fe hyperfine field of a;=0, and 0.048(3) — the y-axes 
in Fig. 31(b) arc already scaled with this factor. 

The spectra of a;=0.145(20) show a single Lorentzian 
absorption line with a temperature independent line 
width in addition to the impurity lines (see above). This 
confirms the //SR results that show the absence of mag- 
netic order of a:=0.145(20) for T>1.6 K. Due to exper- 
imental limitations we did not study the Ce magnetic 
order with Mossbauer spectroscopy because it occurs at 
even lower temperatures (see Sec. , p. 24). 

The low temperature value i?;;/ (12K)=5.6(1) T in the 
undoped compound is reduced to 2.0(1) T for a;=0. 042(2) 
and to 1.0(1) T for x=0.063(2) (see Fig 31). Qualita- 
tively, this is consistent with the reduction of the low 
temperature value of the muon spin precession frequency 
/^(O) shown in Fig. 28. However, is reduced by 75% 
from 28 MHz (a;=0) to 7 Mhz (a:=0.048(3)) but the Fe 
hyperfine field only by 65% from 5.6 T to 2 T if we ne- 
glect the low temperature enhancement. This suggests 



either that the muon site is slightly different in the doped 
materials [47], that the polarization of the Ce moments by 
the Fe sublattice is oriented differently [47], or that the Fe 
magnetic structure is slightly different. 

Below r«10 K, Bhf of a;=0.048(3), and 0.063(2) in- 
creases to ~A T. This increase of the Fe hyperfine field 
does not cause the muon spin precession frequency in 
Fig. 28 to increase. Instead the /iSR transverse relax- 
ation rate At, shown in Fig. 31 increases simultaneously 
with Bhf- This indicates that the additional magnetic 
moment that enhances Bhf is static but disordered. It 
could be due to an increase of the Fe magnetic moment or 
a transferred hyperfine field due to short range Ce order. 



Summary 

We have studied the electronic phase diagram of 
CeFeAs0i_2;Fi: by detailed resistivity, synchrotron X-ray 
diffraction, Mossbauer spectroscopy, magnetic suscepti- 
bility and muon spin relaxation experiments. All phase 
transition temperatures are summarized in Tab. HI. 
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FIG. 32. The electronic phase diagram of CeFeAsO i_ a; F^, 
with a region where magnetic short range order (SRO) and 
superconductivity occur at the same time in the tetragonal 
structure, denoted as "sro+SC+t" At present it is unclear 
whether this is phase separation or microscopic coexistence. 
The close proximity of the transition temperatures Tsro and 
Tc for a;=0. 063(2) indicates that this composition is close to 
the tetra critical point of the phase diagram. All lines are 
guides to the eye. Dotted lines are extrapolations of the phase 
boundaries indicated by solid lines. For a;=0. 150(20) we did 
not study the magnetic properties for T<1.6 K, therefore it 
remains unclear, if the Ce magnetic order survives for this 
composition. 



General results The phase diagram is shown in 
Fig. 32. We find that the structural phase transition 
always precedes (Ts>Tsro) short range magnetic order 
and is preempted by structural fluctuations. The short 
range order develops over a large temperature interval 
and only close to the temperature where nearly 100% 
of the sample volume exhibits magnetic order, we find 
long range magnetic order for T<TN<Tsro- F doping 
suppresses both phases, leading to lower transition tem- 
peratures and smaller order parameters. Seemingly, the 
F doping suppresses the structural transition completely 
before it induces superconductivity. The magnetic or- 
der, however, survives and we find a mixed phase of bulk 
superconductivity and bulk magnetic short range order 
for a;=0. 063(2). For this composition the magnetic order 
never becomes long ranged — whether this is due to the 
competition with superconductivity or simply the same 
short range order that we find for a;=0, and 0.042(2) re- 
mains unclear. Higher doping levels quench the Fe mag- 
netic order and give rise to pure superconductivity. The 
magnetic penetration depth of both superconducting ma- 
terials a::=0. 145(20), and 0.150(20) is consistent with a 
single s-wave gap without nodes. The F doping also 
weakens the Ce magnetic order: As the Fe order becomes 
short ranged for a;=0. 063(2), also does the Ce order and 
for x=0. 145(20) it is suppressed to T^'^sal K — a reduc- 
tion by 75% compared to the undopcd material — and the 



ordered moment is drastically reduced. This suggests 
that the Ce-Ce RKKY interaction is weakened either by 
changes of the Fermi surface, or by the smearing of the 
conduction bands by doping induced disorder. 

Mixed phase That the magnetic order coexisting 
with superconductivity is rendered short ranged for 
x=0. 063(2) is consistent with Ref. 8. We find mag- 
netic order occurs at temperatures slightly lower than 
the critical temperature Tc of superconductivity. This 
allows us to conclude bulk superconductivity close to Tc 
from /iSR. At lower temperatures the /xSR signal is dom- 
inated by magnetic order. However, the temperature in- 
dependent large diamagnetic (macroscopic) susceptibility 
at low temperatures indicates that the superconducting 
volume fraction remains unchanged, i.e., that it remains 
bulk superconductivity also at lower temperatures. Ad- 
ditionally the changes of the lattice parameters below Tc, 
which are related to the pressure dependence of Tc, are 
consistent with those of the two pure superconductors 
studied by X-ray diffraction, which further substantiates 
the evidence for bulk superconductivity. 

The evidence relevant for the question of (nanoscopic) 
phase separation vs. microscopic coexistence is therefore 
similar to the Sm system: 



bulk superconductivity (from /iSR, x(T), 
XRD), 



and 



• bulk short range magnetic order (from /iSR, and 
Mossbauer spectroscopy), 

In sum, above evidence are as ambiguous as for 
the Sm system. We cannot argue with certainty for, 
or against nanoscopic phase separation, or microscopic 
coexistence — the experimental evidence are equally con- 
sistent with both situations: For nanoscopic phase 
separation [7-9], at rsa2 K //SR would show bulk mag- 
netic order, while probes, such as susceptibility measure- 
ments, that are insensitive to the magnetic order, would 
show bulk superconductivity — consistent with our data. 
However, the same evidence would be expected for mi- 
croscopic coexistence. We note that our experimental 
data confirm and extend previously published ^SR and 
magnetic susceptibility measurements [7, 8]. We merely 
argue, that the available evidence are not sufficient to 
decide this matter. 

Relevance for a possible QCP The possible existence 
of a magnetic quantum critical point (QCP) is connected 
to the above discussion. The first phase diagram of this 
system suggested the existence of such a QCP. However, 
Zhao et al. [4] inferred their data from neutron diffrac- 
tion. This technique is, in principle, not sensitive to short 
range order. Taking this into consideration, Zhao's phase 
diagram is corrected by the /LtSR work of Sanna et al. [7] , 
Shiroka et al. [8], and the present work. These findings, 
as discussed above, do not clarify the existence of a mag- 
netic QCP because it remains unclear whether the phase 
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FIG. 33. Intensity maps of X-ray diffraction patterns that illustrate the tetragonal to orthorombic structural transition of 
CeFeAs0i-2,Fa; and its relationship to the magnetic short and long range order, and anomalies in the electrical resistivity. In 
the doped materials the pronounced anomalies of the resistivity disappear and an estimation of Tjv and Ts is not possible by 
the resistivity measurements. 
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TABLE 111. Phase transition temperatures for the structural [Ts obtained by XRD), superconducting (Tc obtained by /iSR), 
Fe magnetic long range order (Tjv obtained by pSR), magnetic short range order {Tsro, Tioo, and AT obtained by fiSK — see 
text). T^"^ is the magnetic ordering transition temperature of the Ce sublattice obtained by /^SR from the fit to asciT), and 
susceptibility measurements (xiT))- AH temperatures are given in K. 
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transition as a function of doping is a first or second order 
transition. 

Electrical resistivity & phase transitions In Fig. 33 
we compare X-ray diffraction, /iSR, and Mossbauer with 
resistivity data. Fig. 33 shows that Ts and T/v of Ce- 
FeAsO can only be estimated from the maximum and 
the inflection point. However, if only the anomalies of 
p{T) were used to determine the critical temperatures, 
T/v would be underestimated and overestimated by 
sa2 K. This is different from SmFcAsO for which the re- 
sistivity maximum is 10 K above the structural transition 
and inflection coincides with T^. Although the cusp of 
p{T) is suppressed in doped samples, we can try to cor- 
relate the resistivity with the structural, and magnetic 
properties. For example, for a;=0. 042(2) the resistivity 
minimum is close to the onset of short range magnetic 
order at T^ro- Similarly, the inflection point of p{T) is 
close to Ts- However, these general correlations do not 
seem to be generic features of p{T) as the structural tran- 
sition is close to the minimum and short range order at 
even lower temperatures for a;=0. 048(3). 

Electrical resistivity & fluctuations Because the max- 
imum of p{T) of a;=0.00 is so close to the structural 
transition, we could speculate about the maximum to 
be caused by the structural fluctuations marked by the 
increase of the FWHM of the (2,2,0)x Bragg peak. How- 
ever, in Sec. , p. 16 wc already learned that for the Sm 
system the increase of the FWHM is not correlated with 
the maximum of p{T). Therefore it seems that Ts^^T^^^^ 
for the Ce system with a;=0.00 is coincidental. The seem- 
ing insensitivity oi p{T) to structural fluctuations and the 
phase transition is even more apparent for x=0. 048(3). 
Here, the structure stabilizes below Ts, yet p(T) in- 
creases. Due to the lack of any systematic correlation 
between p{T) and the FWHM, and the transition tem- 
peratures we have to conclude, as for the Sm system, 
that the structural and magnetic fluctuations studied in 
this work have no simple connection to the temperature 
dependence of the resistivity [43, 98]. 



CONCLUSION 

In Fig. 34 we show a superposition of the two phase 
diagrams. The detailed discussions of the phase diagrams 
can be found on p. 16, and p. 27. 

General results Both parent materials (CeFeAsO and 
SmFeAsO) show the well known structural, Fe magnetic, 
and rare-earth magnetic phase transitions. The F doped 
materials also undergo these phase transitions, but the 
ordering temperatures and order parameters are reduced 
by doping. The reduction of the phase transition temper- 
atures by doping is accompanied by the disappearance of 
the resistivity anomalies. For the studied compounds, 
the disappearance of the resistivity anomalies by doping 
is comparable in both systems. 
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FIG. 34. Superposition of the phase diagrams of 
iiPeAsOi-iFa; with R=Sm, and Ce. Tjv, and Ts seem to 
follow the same trend in both systems under doping, whereas 
Tsro is smaller than Ts in the Ce case and larger than Ts 
in the Sm case. This suggests that Tsro does not reflect an 
intrinsic energy scale but rather the strength/abundance of 
disorder or (micro)strain in the samples. However, the fact 
that short range order occurs for the two samples that show 
a mixed phase of superconductivity and short range mag- 
netic order (1=0.06, and 0.063(2) for _R=Sm, and Ce) sug- 
gests that the clean systems might exhibit a low temperature 
long range magnetic order in the superconducting phase. The 
phase transition temperatures are listed in Tab. I, p. 19, and 
Tab. Ill, p. 29. 



Structural transition The structural transition tem- 
perature and order parameter are reduced by doping. It 
is absent once superconductivity is induced by doping in 
both systems (a;=0.06 for the Sm, and a;=0. 063(2) for 
the Ce system). For the Sm system, this is supported by 
the absence of any additional anomaly of the tempera- 
ture dependence of the speciflc heat. In agreement with 
the X-ray diffraction study of Margadonna et al. [11], our 
experimental evidence disprove the claims of a structural 
transition at T'-^^ITO K in superconducting Sm materials 
by Martinelli et al. [12] We suggest that the FWHM of 
a structural Bragg peak should only be considered as ev- 
idence for a structural phase transition if supported by 
additional experiments. Although the broadness of the 
structural transition (see below) hinders the determina- 
tion of the critical exponents we were able to extract the 
following trend: Doping causes the universality class of 
the structural transition to change from 2D-Ising (un- 
doped) to 3D-Ising (doped close to superconductivity). 
This is similar to the reports of a universality class change 
of the magnetic transition [24]. 

Superconducting gap symmetry We analyzed the tem- 
perature dependence of the in-plane magnetic penetra- 
tion depth of SmFeAsOo.9oFo.io and CeFeAsOi-^Fx 
with a;=0.063(2), 0.145(20), and 0.150(20) and found 
that a single s-wave gap without nodes agrees best with 
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the data. However, within the experimental uncertainty 
a second gap cannot be excluded. In Fig. 35 we plot Tc 
vs. 1/A^^. The here measured superconductors roughly 
follow the linear relationship [99] of other doped i?FeAsO 
superconductors. 

Mixed phase Once doping induces superconductivity 
we do not find any evidence for a structural transition 
(a;=0.06 for the Sm system, and a;=0. 063(2) for the Cc 
system). But we cannot exclude that at slightly lower 
doping levels superconductivity can be found in the or- 
thorhombic phase [11]. Instead we find a mixed phase of 
superconductivity and short range magnetic order. Both 
materials have in common that both phases have nearly 
identical transitions temperatures and that superconduc- 
tivity is bulk directly at the transition but the magnetic 
order only becomes bulk at ~5 K. We analyzed the struc- 
tural data, electrical resistivity, micro- and macroscopic 
magnetic susceptibility, and the magnetic volume frac- 
tion of this mixed phase. Based on our results and the 
results of Refs. 7, 8, and 100 we conclude that, it re- 
mains unclear whether this mixed phase is nanoscopic 
phase separation (two thermodynamic phases) or micro- 
scopic coexistence (one thermodynamic phase). This is 
mainly caused by the difficulty to detect a possible cou- 
pling of the order parameters if both ordering tempera- 
tures are very close to each other. It is much clearer in, 
e.g., K doped BaFe2As2 for which phase separation[55] 
and microscopic coexistence [86] are unambiguously dis- 
tinguishable. 

Relevance for a possible QCP The phase diagram of 
the Ce system by Zhao et al. [4] in principle allows a 
magnetic quantum critical point where Tis!=Tc=Q. This 
feature was already clarified by the ^SR works in Refs. 8 
and 9 and our work confirms these results: both the Sm 
and the Ce system show a region in the phase diagram 
in which magnetic order and superconductivity occur in 
a mixed phase. 

Our study cannot give any direct insight into the possi- 
ble existence of a nematic QCP [91]. The nematic phase is 
mainly marked by electronic anisotropies (see e.g. Ref. 45 
and references therein) in the a-b plane and its study 
should require large enough single crystals. However, 
the nematic order parameter and the orthorhombic dis- 
tortion are directly connected [43]. We showed that the 
low temperature limit of the orthorhombic distortion 5q 
depends linearly on the transition temperature Ts (see 
Fig. 24). This proportionality is the same for the Ce and 
Sm system and is also consistent with published data 
[4, 11]. Margadonna et al. [11] reported two supercon- 
ducting materials that undergo the structural phase. 
of these materials is seemingly doping independent and 
much higher than expected from the linear extrapolation 
oi5o{Ts) (while Ts continues to decrease). This suggests 
that the structural phase transition temperature and or- 
der parameter do not simultaneously vanish; in contrast 
to the simultaneous vanishing of Ts and 5q expected at 



a structural QCP (or the related nematic QCP[91]). 
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FIG. 35. The scaling of Tc with the superfluid density 
A^^ oc Us- CeFeAsO 1- a; Fa: and SmFeAsOi-iFj, supercon- 
ductors follow a linear dependence similar to other reports. 
This is an empirical relationship [99] also found for the hole 
and electron doped cuprate superconductors and other exotic 
materials. Data for other ferropnictides and cuprates have 
been taken from references [5, 6, 65, 67, 101]. 

Broadness of phase transitions Both, the magnetic 
and structural transitions of both systems are signifi- 
cantly broadened. Evidence for the broadening of the 
magnetic transition is the slow increase the magnetic vol- 
ume fraction. The significant broadening of the (2,2,0)t 
tetragonal Bragg peak and subsequent convex tempera- 
ture dependence of the orthorhombic distortion indicate 
that also the structural transition is broadened. The es- 
sential difference between the two systems is that in the 
Sm system the magnetic long range order develops in 
clusters and precedes the structural transition {Tsro>Ts) , 
whereas in the Ce system clusters of magnetic short range 
order occur after the structural transition (Tsro<Ts)- 
Compared to the Ce system, the Sm materials have a up 
to 70% larger FWHM of the (2,2,0)t tetragonal Bragg 
peak. This could mask the onset of the orthorhombic 
order in the small magnetically ordered regions. 

Reports of broad transitions exist, but sometimes are 
not explicitly mentioned in other studies. Sauna, Shi- 
roka, and co-workers [7-9] reported ^SR data and used 
a Gaussian distribution of magnetic transition temper- 
atures with a standard deviation crwS-lO K to quan- 
tify the transition width. This indicates a total width 
of the magnetic transition from onset to saturation of 
Tioo - rsro«4cr«20-40 K (95% of a Gaussian distribu- 
tion lie within ±2(7 of the average). Consistent with 
a broadened structural phase transition, Qureshi et al. 
[50] reported a "precursor" of the orthorhombic struc- 
ture already above the structural phase transition of 
LaFcAsOi_a;Fa; that manifested in an increase of the 
FWHM of the studied Bragg peaks. Considering the 



32 



existing reports of broad structural and magnetic tran- 
sitions, the broadness should be considered as a mani- 
festation of an intrinsic property. The spin fluctuations 
preceding the magnetic phase transition should, via the 
magneto-elastic coupling, play an important role for the 
broadening of the FWHM of the Bragg peak. In addition, 
(random) stress, and disorder may also contribute to the 
broadening of both phase transitions [43, 44, 83, 85]. 

Phenomenological model for the broadened magnetic 
transition The broadening of the Fe magnetic transi- 
tion is characterized by a power law increase of the mag- 
netic volume fraction, a peak of the transverse relaxation 
rate (width of the distribution of ordered moments/local 
magnetic fields) and an unusual temperature dependence 
of the measured average ordered magnetic moment, the 
muon spin precession frequency. We reproduce these fea- 
tures in a simple model for the Fe magnetic transition. 

The model assumes that for T^Tgro (the onset of 
the magnetic order) small, magnetically ordered regions 
form, that have an independent temperature dependence 
of the order parameter. All regions have the same criti- 
cal exponents, but a full ordered moment (at T=0) that 
depends on the local ordering temperature. The only 
adjustable parameters in this model are the critical ex- 
ponents. 

This model fully describes the magnetic ordering of 
SmFeAsO as seen by /zSR. The doped Sm materials show 
an additional contribution to the transverse relaxation 
rate that is most likely caused by enhanced disorder. 
The Ce system does not show muon spin precession for 
TM<T<Tsro but all other features are also present at the 
magnetic transition. This suggests, that the same mech- 
anism is present in the Ce system, albeit with a shorter 
magnetic correlation length. 

Magnetic order of the rare earth Both systems show 
rare earth magnetic order. The most accurate determi- 
nation of the ordering temperature is possible with spe- 
cific heat measurements. The ordering temperature and 
magnetic correlation length is reduced by F doping. /xSR 
indicates that Sm order is present also close to optimal 
doping with a transition temperature that is only reduced 
by 1-2 K compared to the undoped material. For the su- 
perconducting a;=0.06 specimen, the Sm magnetic order 
is short ranged, just as the Fe magnetic order at this dop- 
ing level. The Ce order is affected by doping much like 
Sm but the transition temperature is reduced to K 
for x=0. 145(20). 

This indicates, that the interaction between the 4/ mo- 
ments is weakened by doping. Considering the decrease 
of the lattice constants under doping (see Figs. 3, 22) 
this excludes super exchange interaction as the main in- 
teraction because it sensitively depends on the distance 
between the interacting moments. RKKY interaction 
between the 4/ moments, on the other hand, may be 
weakened by doping, which could be attributed to the 
disappearance of the Fe magnetic order, and changes 



of the Fermi surface topology [96]. Additionally, the 
disorder induced by doping may smear out the con- 
duction bands which effectively leads to a weakening 
of the RKKY interaction. The smearing of the con- 
duction bands was studied for transition metal doped 
BaFe2As2 and CaFe2As2[102]. The effect for F doping of 
i?FeAsO should be weaker, but the additional electronic 
inhomogeneities[103] that develop in these systems may 
play a similar role. 

Rare earth spin fluctuations Earlier reports associ- 
ated the magnetic fluctuations in the Sm system with Fe 
spin fluctuations preceding superconductivity [68]. Our 
detailed analysis of the magnetic fluctuations indicates 
that this is completely coincidental and that the fluctua- 
tions seen by /zSR should be associated with fluctuations 
of the Sm moments due to crystal electric field effects. 
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